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1. Introduction 



One of the main features of analysis on post-critically finite self-similar (pcfss) sets is 
that it is possible to understand the behavior of the Laplacian and its inverse, the Green 
operator, in terms of the self-similar structure of the set. Indeed, a major step in the ap- 
proach to analysis on self-similar fractals via Dirichlet forms was Kigami's proof [8, 10] 
that for a self-similar Dirichlet form the Green kernel can be written explicitly as a series 
in which each term is a rescaling of a single expression via the self-similar structure. In [6] 
this result was extended to show that the resolvent kernel of the Laplacian, meaning the 
kernel of (z - A) -1 , can also be written as a self-similar series for suitable values of z e C. 
Part of the motivation for that work was that it gives a new understanding of functions of 
the Laplacian (such as the heat operator e' A ) by writing them as integrals of the resolvent. 

The purpose of the present work is to establish estimates that permit the above approach 
to be carried out. We study the functions occurring in the series decomposition from [6] 
(see Theorem 3.3 below for this decomposition) and give estimates on their decay. From 
this we determine estimates on the resolvent kernel and on kernels of operators defined as 
integrals of the resolvent kernel. In particular we recover the sharp upper estimates for the 
heat kernel (see Theorem 10.2) that were proved for pcfss sets by Hambly and Kumagai [4] 
by probabilistic methods (see also [1, 13, 3] for earlier results of this type on less general 
classes of sets). It is important to note that the preceding authors were able to prove not 
just upper estimates but also lower bounds for the heat kernel, and therefore were able 
to prove sharpness of their bounds. Our methods permit sharp bounds for the resolvent 
kernel on the positive real axis, but we do not know how to obtain these globally in the 
complex plane or how to obtain lower estimates for the heat kernel from them. Therefore 
in this direction our results are not as strong as those obtained in [4]. However in other 
directions we obtain more information than that known from heat kernel estimates, and we 
hope that our approach will complement the existing probabilistic methods. In particular 
we are able to obtain resolvent bounds on any ray in C other than the negative real axis 
(where the spectrum lies), while standard calculations from heat kernel bounds only give 
these estimates in a half-plane. 

A further consequence of our approach is that we extend (in Theorem 9.7) the decom- 
position from [6] to the case of blowups, which are non-compact sets with local structure 
equivalent to that of the underlying self-similar sets. The blowup of a pcfss set bears the 
same relation to the original set as the real line bears to the unit interval, see [17] for details. 

The structure of the paper is as follows. In Section 3 we recall some basic features 
of analysis on pcfss sets, as well as the main result of [6], which is the decomposition 
of the resolvent as a weighted sum of piecewise eigenfunctions. Section 4 then discusses 
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the natural decomposition of the pcfss set according to a self-similar harmonic structure. 
The next two sections deal with sharp estimates on the real axis. In Section 5 we obtain 
estimates of piecewise eigenfunctions for which the eigenvalues are real. To show these we 
decompose into pieces in the manner of Section 4, prove that each piece has one normal 
derivative that is larger than the others by a factor, and use this to show that the only way 
to glue pieces together smoothly is by requiring the size to decay by this factor each time 
we cross a cell of the decomposition. This implies that piecewise eigenfunctions have sub- 
Gaussian decay. From these estimates it is routine, though somewhat long and technical, 
to estimate the kernel of (A - A) for A on the positive real axis and to prove that the 
decomposition of [6] is valid on blowups in this setting. The former is in Section 6 and the 
latter in Section 7. 

In Section 8 we switch gears and prove some complex analytic estimates related to the 
Phragmen-Lindelof theorem. Our goal is a method for obtaining resolvent estimates on 
sectors in the complex plane from our decay estimates on the real axis. This method is not 
restricted to the setting of pcfss sets, and may be of interest for proving resolvent estimates 
in more general settings, such as metric measure spaces. We find that combining decay 
estimates on the real axis with some weak estimate on a sector away from the spectrum 
proves decay on any ray in the sector. In Section 9 we show that the required weak esti- 
mates can be proved on pcfss sets by a generalization of some arguments from Chapter 4 
of [9]. As a result we obtain our main result, Theorem 9.6, giving decay estimates for 
piecewise eigenfunctions and the resolvent kernel away from the negative real axis in C. 
Section 10 contains some examples, including a proof of the upper bounds for the heat 
kernel. 

2. Acknowledgements 

This paper relies heavily on its precursor [6]; the author thanks his co-authors on that 
paper, especially Bob Strichartz, whose suggestion that there might be an approach to heat 
kernel estimates via Kigami's theory and the results of [6] inspired the present work. 

3. PCFSS Fractals: Energy, Laplacian, Metric, Eigenvalue Counting and Resolvent 

KERNEL 

We work on a pcfss fractal X corresponding to an iterated function system [F\, . . .Fj}. 
In full generality this could be a compact metrizable topological space X equipped with 
continuous injective self-maps Fj so that there is a continuous surjection from the space of 
infinite words over { 1 , . . . , 7} N to X, as in Chapter 1 of [9] . The reader may prefer, however, 
to think of the more intuitive situation in which the Fj are Lipschitz contractions on a finite 
dimensional Euclidean space and X is their fixed point in the sense of Hutchinson [5]. 
In either case, the boundary of X is the post-critical set Vq, which is a finite set with the 
property that Fj(X) n F k (X) c F/V ) D F k (V ) for k * j. We assume that X \ V is 
connected. The best known example of a set of this type is the Sierpinski Gasket, for 
which the full analytic theory we will use is expounded in [19]. Note that the Sierpinski 
Carpet is not of this type, as it is not post-critically finite. 

If w is a finite word on the letters { 1 , . . . , J} then |w| is its length. The set V m is 
U| M ,|- m F M ,(Vo) and V* = U™ =0 V m . Points in V* \ Vq are called junction points. For a word 
w = w\ ■ ■ ■ w m we write F w = F m o F Wt o ■ ■ ■ o F Wm , The set of finite words is W t , and 
the set of infinite words is S. If cr — cr\(j2 ■ ■ ■ is an infinite word we use [<x] m to denote 
the subword cr\ ■ ■ ■ cr m . For a finite word w we use the notation wW t (respectively wL) for 
finite (respectively infinite) words that begin with w. 
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We assume that we have a regular self-similar harmonic structure on X which provides 
an irreducible self-similar Dirichlet form E with domain dom(£) and scalings < r j < 1 
so that 

j 

E(u) = Y u rfE{uoF j ). 
i 

Details of the construction and properties of such forms may be found in [9]. We also fix a 
probability measure p satisfying 

j 

(3.1) f i(A) = Y J W( F ] 1 ( A »- 

i 

We will require that this measure is related to the Dirichlet form by pj = r? for the unique 
S so that Yjj 'j — 1 ■ This measure is known to be the natural one for many aspects of the 
analysis on X; it arises in studying the analogue of Weyl-type asymptotics of Laplacian 
eigenvalues [12], and in determining heat kernel estimates [9, 13, 3, 4]. It is the Hausdorff 
measure (with Hausdorff dimension S ) for the resistance metric, which is the metric in 
which the distance between points x and y is R(x, y) given by 

flOt,);) -1 = min{£ , (w) : u(x) = 0, u(y) = 1}. 

Remark 1. In most of the theory of analysis on fractals one can take any Bernoulli measure 
defined by (3.1) for some < pj < 1 with Yij^j - 1- In particular, this is all that is required 
for the construction of the Laplacian resolvent in [6]. The author would be curious to 
know whether it suffices for Theorem 9.7 or even Theorem 7.3 regarding the resolvent on 
blowups. 

The Dirichlet form and measure give rise to a weak Laplacian A by defining u e dom(A) 
with Am = / if there is a continuous / such that 



E(u,v) = - J fvdp 

for all v e domo(£'), the functions in dom(£) that vanish on Vq. The Laplacian is self- 
adjoint and has compact resolvent, with negative eigenvalues -Aj of finite multiplicity that 
accumulate only at -oo. The asymptotic distribution of the eigenvalues was determined by 
Kigami and Lapidus [12], who proved a more general version of the following result. 

Proposition 3.1 ([12]). Let N(x) = #{Aj : \Aj\ < x). Then 

We also recall from [9] that there is an explicit formula for a continuous Green kernel 
that is positive on the interior of X, zero on Vq, and inverts -A with Dirichlet boundary 
conditions. Kigami's construction was generalized in [6], from which we will need the 
following results. 

Proposition 3.2. If p e Vq and z e C is such that none of the values r w /j w z is a Dirichlet 
eigenvalue of A then there is a function r\$ with measure -valued Laplacian satisfying 

f(z/ - A)j#> =0 on X \Vo 
rjf = S pq for qeVo 

where 5 pq is the Kronecker delta. We call this a piecewise z-eigenfunction. 
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Under the same assumptions on z we define functions that form a natural basis for the 
piecewise z-eigenfunctions on 1 -cells. For p e V\ \ Vo, let 

(3-2) ^-X'r^'v 

where by convention the sum ranges only over those terms that are well-defined, in this case 
those j so p e Fj(Vo). Then tfrp has a measure-valued Laplacian, with (zl - A)^/f = 
on all 1 -cells and Dirac masses at the points of V\ \ Vq. The main result of [6] is that 
for suitable values of z the resolvent G^ z \x,y) of the Laplacian with Dirichlet boundary 
conditions may be written in terms of these piecewise eigenfunctions. 

Proposition 3.3. For each z e C such that none of the values r w p w z is a Dirichlet eigen- 
value of A, let 

(3.3) G (z \x,y) = J] r w ^" z \F- l x,F- l y) 

weW, 

in which 

(3.4) ¥ z \x,y) = ^ G^\x)^f(y) 

where Gpl is the inverse of the symmetric matrix Bpl with entries 

(3.5) Bf q = ^ d^^ig). 

j:Fj(Vo)sq 

Then G^(x,y) is continuous and vanishes on Vo, and it inverts (zl — A) in that if f 6 L l (fi) 
then 

(Zl-A) f G (z \x,y)f(y)dp(y)=f(x). 
Jx 

In what follows, we will suppose that z = A > and study the decay of r] w and its 
dependence on A, so as to estimate the decay of the terms in G (1) . This will then be used to 
study the dependence of rj^ and G (:) for those z € C that are not on the negative real axis. 
As usual these estimates will include constants depending on the structure of the fractal, 
for which reason we define a shorthand notation as follows. 

Notation 3.4. We write a < b or b > a if there is a constant c > that depends only on the 
fractal, the harmonic structure or the measure, and for which a < cb. If a < b and b < a 
then we write a b. The implicit constant, as well as constants explicitly named, may 
vary from line to line in a computation. 

4. Resistance Partitions and Chemical Paths 

The short-time behavior of diffusion on fractal sets may be analyzed by partitioning the 
fractal such that all pieces have a prescribed resistance (up to a constant factor) and study- 
ing the lengths of minimal paths in this decomposition (see, for example, [4]), which are 
often called chemical paths. We should therefore expect that the behavior of the resolvent 
G^ A \x,y) for large A > may be determined in the same manner. In this section we record 
some known estimates for chemical paths that will be useful later. 

Definition 4.1. A partition of X is a finite set 0cff, with the properties 

(1) If 0,0' e® then n 6"L + only if = 0', 

(2) U ee0 6E = 2. 
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We will refer to the cells F 9 (X), 9 e 0, as the cells of ©. 
The partitions we use are 

®k = {Ol ■ ■ ■ dm ■ r 9l ■ ■ ■ r Bm < e~ k < r Bl ■ ■ ■ r 6m ^ }, 
and we write V(®a) for the junction points corresponding to words in ©a 

V(® k ) = (J Fg(V Q ). 

We will frequently use that when 9 e ©a, 

(4.1) r e ^e-\ Ho^e- ks , We - e -* (s+1) . 

There is a graph structure corresponding to the partition ©a. Let be the graph with 
vertices V(©a) and edges between every x and y for which there is 9 e ©a such that 
x,y 6 Fg(Vo). Points joined by an edge are said to be neighbors. As a consequence of 
the definition of ©a we have (see Lemma 3.2 of [4] for a proof) 

(4.2) R(x,y) - e~ k for any x,y e V(©a) that are neighbors in rV 

A path joining x and y in Fa is a finite sequence fe}[_ of vertices such that adjacent 
vertices are connected by an edge and xq = x, Xi = y. The length of the path is / and the 
graph distance dk(x,y), also called the chemical distance, is the length of the shortest path 
joining x and y in r*. 

It is generally difficult to obtain good estimates relating the graph distance dk(x, y) and 
the resistance R(x,y); the interested reader is directed to [11] for a detailed analysis of 
this question and its connection to heat kernel estimates. We will satisfy ourselves with 
some elementary but crude results. For example, it is well known that R{x,y) does not 
exceed the resistance along a path from x to y, so that for x,y € V(©a) we must have 
R(x,y) < e~ k dk{x,y). In particular if x, y e Vo then dk{x, y) > e k . Conversely, in Lemma 3.3 
of [4] it is shown that for any x,y € X, the bound dk(x,y) < g( s+ W 2 holds. We will need 
these estimates in a form that compares dk+k'(x,y) to dk(x,y). 



Lemma 4.2. 



e k' < dk+k'{x,y) (s+iyt/2 
d k (x,y) 



Proof. Observe that the partition by words 

{9" = 99* :9e& k ,9'e & t ] 

has rg» < e~ (k+k \ so is a subpartition of ®k+k> and consequently has longer paths. Since a 
path in this partition describes a path F^ in which any pair of vertices is joined by a path 
no longer than the maximal path between points of Vo in F^ , which we know is at most 
e (s+\)k'/2^ we conclude that d k+k >(x,y) < e (s+1)k ' ,2 d k (x,y). 

For the lower bound we note that there is a constant c depending only on the harmonic 
structure and such that the partition by words 

{9" = 99' :9<E@ k ,0'& ®k>-c] 

contains ®k+k> and therefore has shorter paths. Again each such path restricts to a path in 
Fa, and now every pair of Fa vertices is separated by a path of length at least the minimal 
distance between points of Vo in F^-c, which is > e > e k ' . a 
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For the special classes of nested and affine nested fractals there are quite precise results 
about chemical distances in [13, 3]. In particular the authors construct a geodesic metric 
on affine nested fractals that is comparable to a power of the resistance (see Proposition 3.6 
and Remark 3.7(2) of [3]), from which the following is easily obtained. 

Proposition 4.3 ([3]). For an affine nested fractal there is y € [g+T' l\ wn ^ cn ma y be 
obtained by solving an explicit optimization problem, such that for all sufficiently large k, 

d k (x,y)^e ky R(x, y y. 

Sketch of Proof. Knowing that R(x, y) 7 is comparable to a geodesic metric we see that 
edges of have geodesic length like e~ Ay , so the geodesic distance from x to y is com- 
parable to d k (x,y)e~ ky if k is large. This must be comparable to R(x,y) r , so we have 
d k (x, y) — e ky R(x, y) y . The bounds on y are evident from the preceding discussion. □ 

The proof in [3] is a generalization of an earlier argument, Proposition 3.5 of [13], for 
the case of nested fractals. We do not wish to give the definitions of nested or affine nested 
fractals here, but we recall that they are subsets of R", are generated by iterated function 
systems consisting of Euclidean similarities that have a high degree of symmetry. Full 
details may be found in [13, 3] and the references therein. 

Definition 4.4. For later use we extend the definition of d k (x, y) from to all of X in 
the obvious fashion. For x in the interior of a cell F e (X) of ® k and y e F^ let d k (x,y) = 
min{d k (z,y) : z e Fq(Vq)}. For x in the interior of F e {X) and y in the interior of F ff {X) set 
dk(x,y) = min{c4(z, z') : z e Fg(Vo), z' e F s ,(Vq)}. When we refer to a d k geodesic between 
points x and y that are not in 0^ we mean a geodesic joining the cells containing them. 

5. Estimates for piecewise eigenfunctions with positive eigenvalue 

In this section we develop decay estimates for the piecewise /l-eigenfunctions tK in the 
case where A is a positive real number. We summarize the results in the following theorem. 

Theorem 5.1. There is k(A) with e kiX) (1 + / i) 1 / (s+1) and a constant c > depending only 
on the fractal and harmonic structure such that for A e (0, oo) and p 6 Vq, each piecewise 
A-eigenfunction rffi satisfies the following bounds. For all q € Vq, q ± p, and all x e X 

(5.1) exp(-cd k (X)(p, x)) < rfpHx) < exp(-d k{A) (p, x)), 

(5.2) d nn f{p) - (1 + A) l « s+l \ 

(5.3) -(1 + ^) 1/(S+1) exp(-d kw (p, q)) < d n ^\q) < -(1 + ^) 1/(S+1) exp(-cd kw (p, q)), 

except that the lower bound of (5.1) is not valid on cells F e (X), 8 € such that F e (Vo) 
contains a point of Vq \ {p}. 

Proof. The number k(A) is introduced in Definition 5.12, where it is shown it is comparable 
to (1 + A) x ^ s+l \ The upper bound of (5.1) for sufficiently large A is Corollary 5.15, and 
continuity ensures we may take a suitably large constant multiple to make it true for all 
A > 0. Both (5.2) and the lower bound of (5.3) are in Lemma 5.18, while the lower bound 
of (5. 1) is Corollary 5.21 and the upper bound of (5.3) is Corollary 5.22 □ 

Remark 2. We note in passing that the estimate (5.1) is sufficient to complete Strichartz's 
proof of the smoothness of finite-energy harmonic functions on products of fractals [18, 
Theorem 1 1 .4] . Strichartz proves this modulo an assumption on the behavior of the normal 
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derivatives of the heat kernel [18, Equation (8.3)], but the key estimate in his proof is ac- 
tually the decay estimate [18, Equation (8.25)], which is a consequence of the above in the 
case of nested fractals by Proposition 4.3. More generally, for products based on fractals 
and harmonic structures of the type discussed in this paper, the estimate (5.1) is sufficient 
to imply the convergence in [18, Equation (8.19)], and thus the argument proving [18, The- 
orem 11.4]. The author hopes that is a first step toward proving hypoellipticity of solutions 
of elliptic PDE on products of pcfss fractal sets ( [18], Section 1 1 and [16] Section 8). 

The remainder of the section is devoted to the proof of this theorem. As the working 
is at times technical it may help the reader to have a concrete example in mind. The most 
elementary choice is to let X be the unit interval with Lebesgue measure and the maps Fj be 
the contractions onto the left and right halves of X. In this case we obtain the usual Dirichlet 
energy, the Laplacian is the second derivative, and the normal derivative is the outward- 
directed first derivative at the endpoints [19]. Eigenfunctions and piecewise eigenfunctions 
are exponentials, and a quick computation shows that rf* = sinh ~\/~Ax/ sinh VX Our goal 
is to show that the basic features of this function, as listed in the theorem, are also true for 
the functions rfy on any pcfss set. 

An observation we use throughout this section is (5.4), which says that we may decom- 
pose these functions into pieces corresponding to cells of a partition ®k and obtain a linear 
combination of pieces; moreover if ©^ is suitably chosen then each piece looks almost like 
a copy of a fixed function. In the case of the unit interval the pieces are literally translates 
and reflections of a fixed hyperbolic sine, and the decomposition (5.4) can be obtained by 
iteration of the double angle formula (see Section 2 of [6]). The main idea of this section 
is that in the decomposition each piece has a larger normal derivative near its peak than at 
the other boundary points, and in order for such pieces to join together in a smooth fashion 
it is necessary that the peak of each is smaller than that of its neighbor by a constant factor. 
Thus the exponential decay may be derived from the shape of a collection of basic func- 
tions and counting the cells in paths on In essence, our estimates on the relative sizes 
of the normal derivatives come from the fact that rfy is subharmonic, but there is some un- 
avoidable technical work to obtain the correct quantative bounds. It should be emphasized 
that most of the work in this section depends on the fact that A is positive. 

Assumption 5.2. In this section we require that A > 0. 

We will extensively use the decomposition of rfy as a piecewise /l-eigenfunction on the 
cells of a partition ©^ as defined in Section 4 above. Using the Laplacian scaling 

A(//r - o f- 1 ) = {r w ^r\^ A) ) o f- 1 = a,;;*' ^ o f- 1 

we see that rffi may be written as a linear combination of the functions 7/ ( v rf,/if,/1) for e ©^ 
and y e V(Qjt)- More precisely, 

(5-4) rjf\x) = J] ti*™ fas'* ° F e 1 

0e@ t qev 

Before proceeding we require several preparatory lemmas about the function rf}' . The 
first is a maximum principle for smooth subharmonic functions. It is well known but does 
not seem to appear in the literature, except for that part which is in Proposition 2.11 of 
[15]. 

Lemma 5.3. Suppose u e dom(A) and Au > 0. If u attains its global maximum at an 
interior point then u is constant. Moreover Au < at any local maximum point. 
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Proof. The proof uses the fact that on any cell F H (X) there is a Green kernel g w that is 
non-negative on F H ,(X) and strictly positive on F W (X) \ F w (Vq), and such that 



where h w (x) is the harmonic function on F W (X) with h w (x) = u(x) for all x e F w (Vo). 

Let h m be the piecewise harmonic function at scale m with h m (x) = u(x) for all x e V m . 
Since Am > we find from (5.5) that u < h m . This implies that d„u(x) > d n h m (x) for all 
x e V m . However the sum of the normal derivatives of u at any point of V m \ Vq must vanish 
because u e dom(A), so the sum of the normal derivatives of h m must be non-positive. This 
gives that the ;7i-scale graph Laplacian of /?,„ is non-negative, so if /?,„ achieves its maximum 
at a point x € V m \ Vq then this maximum is also attained at all neighbors of x in the ra-scale 
graph. Connectivity then implies h m is constant on V m . 

If u attains its global maximum at a point x e X \ Vq then u < h„, and h m harmonic 
implies that h m achieves the same value at a point on the boundary of any cell containing 
x. For all sufficiently large m, this point (which could be x) is not in Vo, so h m attains 
its maximum value u(x) at a point of V m \ Vq, and is therefore constant by our previous 
reasoning. Applying this for all large m implies u = u(x) on the dense set V*, and since u 
is continuous it must be constant. 

Suppose in order to obtain a contradiction that u has a local maximum at x e X \ Vq and 
Am(jc) > 0. The easy case is when x e V*, because it is then easily seen that h„, has a local 
maximum at x for all sufficiently large m. This contradicts the above reasoning showing 
the graph Laplacian of h m to be non-negative if there is a neighborhood on which Am > 0. 
The alternative is that x i V t . Then there is an infinite word <x e £ so n m F[o-] m (X) = \x] 
and these sets form a neighborhood base of x. From (5.5) we see that h[cr] m (x) > u(x) for all 
m large enough that Am > on F[a-] m {X). The maximum principle for harmonic functions 
then gives 



so that every neighborhood of x contains a point at which u exceeds u(x), in contradiction 



Definition 5.4. Let £ p be the function harmonic on X with (, p {q) = for q e Vo, q + p and 

Up) = i- 

Lemma 5.5. < rf*> < C P 

Proof. If the first inequality fails then the fact that rfp^ > on Vo implies that tyy has a 
strictly negative minimum at some interior point x. However by Lemma 5.3 the Laplacian 
must be non-negative at a minimum point, in contradiction to Arfp (x) = Atj p (x) < 0. Hav- 
ing established the first inequality, it follows that r/y is subharmonic and thus is bounded 
above by the harmonic function with the same boundary values, which is precisely £ p . □ 

Corollary 5.6. On any connected open set in X the maximum of rfy is attained at the 
boundary. 

Proof. We have Ar] p = Arfp > and so may apply Lemma 5.3. □ 
Lemma 5.7. If A > A' then iff < iff\ 



(5.5) 




ma\{u(y) : y e F W (V )} = ma\{h [cr] Jy) : y e F W (V )} > u(x), 



to u(x) being a local maximum. 



□ 
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Proof. Supposing the contrary we find that rf p ' - rf p has a negative local minimum at a 
point x. However using that rf p 1 > and the hypothesis, 

A(^ n - i/^X*) = (A'i^ - AnfXx) < Xrfp - !#>)(*) < 
which cannot occur at a minimum point by Lemma 5.3. □ 



Corollary 5.8. For any p, q € Vq with p + q, the values d„T] p (p) are non-negative and in- 
creasing in A, while the values d„T) p (q) are non-positive and increasing in A. In particular 
for p + q the values d n Ty p (q) are bounded below by min Mg v Q d„( p (q) < 0. 

Proof. Positivity of d n T] p (p) is evident from Lemma 5.5 because rj p < Cp < 1 on X \ \p). 
Then the fact that rf^ < rf p ^ when A > A' from Lemma 5.7 implies d n rf p \p) > d n rf p \p). 
Similarly d„t] p (q) is negative because rf p > = rf p \q) from Lemma 5.5, and rf p < rf p * 
implies d„T] p (q) > d„rf p \q). The lower bound comes from Lemma 5.5, because rf^ < £ p 
for all A > implies d„rf p \q) > d„( p (q) for all p, q e Vq. □ 



With these basic observations in hand we look more closely at the decomposition (5.4). 
Lemma 5.9. For all p e Vo, 

Jx 

Proof. For A < 1 the result is clear from the positivity of rf p and the monotonicity shown 
in Lemma 5.7; in fact we have rather than an inequality. For A > 1, fix k such that 
e -k(S+r>^ e (e~ l , 1], so that for 6 e 0^ we have rgpgA 1. The monotonicity of Lemma 5.7 
then implies 

f rf rmX> dp. * 1. 
Jx 

All terms in the decomposition (5.4) are positive by Lemma 5.5, so writing 6 P for a word 
in ©a such that Fg p (X) 3 p, we have 

Jx Jf Si ,{X) 

C (.?„ Mo X) 
JX 

However fig ^SKS+D for g e &k and 

our choice of k, so the proof is complete. □ 



Lemma 5.10. For any word w e W*, 

A 



xeF„(y ) {(r w p H ) 1 + A)'" 
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Proof. Applying the Gauss-Green formula to Tfq™ 1 ^" o F n l and the constant function 1 
yields 

= A f rj^oF-Ufi 

Jf w (x) 

Jx 

> Aju w (l + r H p w A)^ 
A 

((r^)- 1 +A)^ 

where we used Lemma 5.9 and = □ 

On the first reading of the following lemma one should think of the case 0' = 0^. It will 
later be used for 0' = 0^. \ {9 : Fg(X) 3 p) and more complicated sets. This lemma is the 
main argument in this section of the paper, in that it uses smoothness of the join between 
pieces of the decomposition (5.4) to show that rfp must decay. 

Lemma 5.11. For 0' c ©t and Y = Ug e @,Fg(X) we have 

£ rjf(y) < c(e k ^A- 1 + £ <>(z). 

yeV(0')\dY zeOY 

Proof. At any point x e V(©a) the sum of the normal derivatives of rfp over the cells 
meeting at x must be zero because rfp e dom(A). If we sum this cancelation over all 
V(©a) points that are interior to Y we may use (5.4) and rearrange to obtain 

(5.6) J] Z ffwM Z d "^ A> ° F o 1 )^) = 0- 

fie©' qeV xeF e (Vo)\dY 

We estimate the innermost sum in (5.6) using Lemma 5.10. Suppose first that 9 and q 
are such that Fgq is an interior point of Y. Then Corollary 5.8 tells us that the terms 
d„{rfq 6 o F g l )(x) for x e dY are negative. Thus the lower bound of Lemma 5.10 is still 
valid with these points removed. Substituting into (5.6) we find 



(5.7) -J] J] 4 A) (^(M Z 9,Mr A) °Fe 1 )^) 

9e@' {qeV :Fi(q)edY] xeF e (y )\dY 

= Z Z ^{Feiq)){ Yu o^Xx)) 



fteQ' («eV :F e ( 9 )«ai'| xeF 8 (V )\57 



* Z Z 



fle© 1 {qeVo:F e (q)tdY} ((rgflg) 1 + X) s 

A 



Z 



Now in the terms 

d n (rf^ x) o F e l ){x) = r^d^^Cx)) 

on the left in (5.7), the points q and Fg l (x) cannot coincide, because x t dY and Fgq e 
dY. Again appealing to Corollary 5.8 we see that such normal derivatives are negative, 
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increasing in A, and bounded below by -c = mm p , ei£ v d„£ p (q), which depends only on the 
harmonic structure of X. Putting this and the estimate r e e~ k into (5.7) gives a positive 
constant C depending on -c and the degree of vertices so that 

C e *2>«( Z )>-£ Z rjfiFdq)) J] (-«*) 

* - Z 7»(y). 

( ce *6S+l) ygy^^y 

We also use that (ce k(s+1) + A)^> < e ks + A^r. □ 

We may obtain decay estimates of 77^ by iterative use of Lemma 5. 1 1 for an appropriate 
choice of k and sets 0'. 

Definition 5.12. Given A > we let k(A) be the larger of and the greatest integer such 
that 

k(A) < j±-j\og A -\ogC -2 
where C is the constant in Lemma 5.1 1. Note that e kw (1 +/l) 1/(s+1) . 

The following result is an immediate consequence of the definition of k(A). 
Corollary 5.13. If A is large enough that k( A) > 1, then for &' c &t(A) 

Z ^^-Yrifiz). 

yeV{&)\8¥ zedY 

Lemma 5.14. For A as in Corollary 5.13 and each p e Vo, let X^ip, - [x e V(®k(A)) '• 
dk(,i)(p, x) > i). Then 

Z T^OO^- 1 

Proof. We induct over i. The base case i — is simply the fact that rjy(p) - 1. Observe 
that points y e dX^)(p, i) for / > 2 satisfy d w (p,y) > i — 1 or are in Vo \ If). The latter 
may be ignored because rjf 1 i 
Lemma 5.13 gives the result. 



may be ignored because rf p is zero on Vq \ {p}, so substituting the inductive estimate into 



If we restrict ify to a cell, Corollary 5.6 implies the maximum is at the boundary. If a 
point x € X has d^ip, x) = i then the boundary points of the cell of ®k(X) that contains x 
are in X^ip, i) Thus Lemma 5. 14 implies 

Corollary 5.15. For any x e X and A as in Corollary 5.13, 

rfp\x) < exp[-d k{A) (p,x)j. 

Lemma 5.14 also allows us to show that Lemma 5.9 gives the correct value for the 
integral as A — » co. 



Corollary 5.16. For all p e V > 



rffdn^il+AT 
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Proof. For A < 1 the result was observed in Lemma 5.9. If A > 1 then for 8 e ®k(A) we 
reason as in the proof of Lemma 5.9 to find 



T] q dp — p.g — A 



iFgQC, 

Integrating the decomposition (5.4) then yields 



X 



9e&(k(A)) qeV a 

= r ^ Z Z 



< /I V e ' 



where the last inequality is from Lemma 5.14. The reverse inequality is Lemma 5.9. □ 

Using the result of Corollary 5. 16 in the proof of Lemma 5.10 improves it as well. 
Corollary 5.17. For any word w € W„, 



From here it is not difficult to obtain estimates of the normal derivatives of t]^ 1 at points 

of Vq. Recall from Section 4 that dk(p, q) is > e k , so exp(— dk(X)(p, <?)) tends to be small 

when A is large. The following lemma therefore tells us that the normal derivatives d n rfp\q) 
are much smaller when q + p than when q = p. 

Lemma 5.18. For p,q e Vq, 

d n ^\p)^{l+A)^, 
-(1 + X)T* exp(-d m (p,q)) < d n rjf(q) < 0. 

Proof. Suppose q + p. Let Q 1 , . . . ,6" be those words from ®k(A) for which F g ,(Vo) 3 q. 
Using (5.4) and the lower bound for the normal derivatives d n (j] r / >i ' fX )(y) for x,y e Vq 
with y + x from Corollary 5.8, as well as r~} e k( ' l> ^ (1 + A) 1/(s+r> , 

> d n rjf\q) = £ Z nfiF^dn^^ o FfXq) 
(=1 .teVo 

= Z Z nf{F, X )r-Mn { ^\F-M) 



>(1+A)^( max -d„Uyj)S] Y.tff>(.F#x) 
>-C(l+X>™ J] rjf(z) 



zeX klA) (g,\) 

> -(1 + exp(l - d m (p,q)) 

where in the final step we applied Lemma 5.14 and used that points in z e Xk(A)(q, 1) have 
dk(,i)(p, x) > dk(A)(p, q) - L This gives the desired estimate for q + p. 
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As a special case of the above estimate we have > d n rf p \q) > —A.***, Combining this 
with Corollary 5.17 for the empty word w gives the desired result for d n T] p (p) when A > 1 . 
For A < 1 the conclusion is clear from Corollary 5.8, because d n j] p (p) is increasing, so is 
bounded below by d„( p (p) and above by d n rf p \p). □ 

To complete our picture of the behavior of r/ (A) we need a lower estimate on its decay 
and a corresponding upper estimate for the normal derivative d„r[p\q) when p + q. After a 
preliminary lemma, these may be obtained by somewhat simpler reasoning than that used 
earlier. 

Lemma 5.19. 

(1) The function rf p is non-zero all points in V* \ Vq, 

(2) For p + qwe have d n rfp{q) < 0, 

(3) For w € ©Ay) and p + qwe have d n t] ( p ntl ' %X) {q) < -1. 

Proof. The first step is to prove a weaker version of the second statement. From the Gauss- 
Green formula we have 

drflf iq) - 8 n ( p (q) = d nn f{q) - d n { q (p) 

= 2 d n Tjf(x){ q (x) - i p l \x)d n ^ q (x) 

xeVo 

Jx 

= A f r,f( q d» 
Jx 

so that d„rfp\q) — > d n ( p (q) as A I 0. With p + q we have d„( p (q) < 0, hence we may find 
some c < and A such that d n rf p (q) < c for all A < A. 

Now at x e V*\ Vo the fact that tj^ e dom(A) requires that the normal derivatives sum 
to zero. Take k so large that x e V(&k) and for all 9 € ©a we have rgpgAo < A. Suppose 
that rfp\x) = 0. Using (5.4) we can write this sum of normal derivatives as 

= Z X ^ w) d »(iT" A) ° F e 1 w) 

l9€& t :F e (y )Bx) p'eVo 

= Z Z nf{F 9 p')r?{d^\F- e \x)) 

{OeQf.FdV^ip'eVo-.Fnip'^x} 

(5-8) <c J] J] r e l ^\F eP ') 

{0e® k :Fe(y o )sx} lp'EV :F e (p')*x) 

where the first step uses the scaling of the normal derivative and the fact that rf p (x) = 0, 
while the second uses that (d„ri p T liA )(F g ~ l (x)) < c because rgjjgA < A and p' + F g l (x). 
However c < and all r^rj^iFgp') > 0, so the only way (5.8) can be true is if these 
values are all zero, meaning that rf p vanishes at the boundary points Fg(Vo) of each cell 
meeting at x. Repeating the argument inductively we see after finitely many steps that rf p 
must vanish at all points in V(0a), which is impossible because p e V(®k) and rf p \p) = 1. 
This proves the first statement of the lemma. 
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The second assertion of the lemma now follows fairly easily from the first. Using the 
partition ®k as above and the decomposition (5.8) for the normal derivative at x = q we see 

d n rj^\q)^ Z 4HF a p') re \d n ^r A) ){F g \ q )) 

\ee® t :F„(V )3q){p>eV a :F (p>)±q\ 

(5.9) <c J] ^ r^fiFep') 

{SE0 t :F 9 (Vo)3?) P'eVo 

and we have already shown that the values r g x rfp\Fgp') > when Fgp' £ Vo- 

Finally the third statement follows from the second because d„rfp\q) is continuous, 
non-positive and increasing in A, and r w fi w A is bounded above when w e ®k(x>, with all of 
these depending only on the fractal and harmonic structure. □ 

Lemma 5.20. There is a constant a > such that for any point x e V(®k(A)) \ Vo> 

r]f{x) >amax{77«00 : 30 e 0, with F 9 (V ) d {x,y}). 

These points y are the neighbors of x in V(®k( A ))- 

Proof. Smoothness of rj^ requires that the first step of (5.8) holds, where we take k = k(A). 
As we are not assuming rf p '(x) = 0we obtain instead of the second step of (5.8) 

rifix) J] r g l d n {if r F r*){F g \x)) 

[Se0 w :F e (V o )3x) 

= " Z Z ^{FeiqWd^r^Fe 1 ^)) 

{6E® m :F e (V a )3x\ {qeV :F 6 (qUx\ 

* Z Z 

where we used that -d„{rlq mX) ){F-\x)) > 1 from the third part of Lemma 5.19. On 
the left of this inequality each of the normal derivatives dn^^'^^Fg 1 (x)) is bounded 

Fy (X) 

above by a constant depending only on the harmonic structure, as shown in the first part of 
Lemma 5.18. Since the values r g l are comparable on both sides of the equation, the result 
follows. □ 

Corollary 5.21. Let X^ip be the subset of X obtained by deleting those cells Fg(X), 
8 e ®k(,i) that intersect Vo at points other than p. There is c > such that for all x € X^),p 

rfp\x) > exp(-cd k(A) (p, x)) 

Proof. If x € Vi®^) \ Vo and d^x){p,x) = i then there is some y e V(®k(A)) with 
dk(A)(x,y) = i — 1. By the previous result, ?] ( p\x) > arfpfy) and by induction rf p X) {x) > 
a'jf p A \p) = a'. The result for points of V(®k(X)) \ Vo follows by setting c = - log a, and we 
note that the estimate Tj ( p\x) < C exp(-dk(A)(p, x)) from Corollary 5.15 ensures c > 0. To 
obtain the bound for a general point x 6 X^AXp, let & £ ®k(A) be such that Fg(X) 3 x. The 
value dk(A)(p, x) is the distance from p to the nearest point of Fg(Vo); all other points of this 
form are at most dk( A ){p, x) + l from p, and none is in Vo \ {/?}, so the result for V(®k(A)) \ Vo 
applies to them. It therefore suffices to know that the restriction of if^ to Fg(X) is bounded 
below by a multiple of its boundary values. Observe that this function is a linear combi- 
nation of the functions rf q r " ,A '' X) o F g 1 with the boundary values as coefficients. Hence it is 
enough to know that each r/f' 1 ^ '(y) has a positive lower bound on those cells Fj(X) that 
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do not contain a point of Vo \ \q). This follows from the third part of Lemma 5.19, conti- 
nuity of t] p and the fact that the values rgpgA lie in a bounded interval, and we see that the 
constant depends only on the harmonic structure of the fractal. □ 

Corollary 5.22. For the constant c of Corollary 5.21 

dnrfp°\q) < -d + A) l/{S+1) exp(-cd k{A) (p,x)) 

Proof. If we rewrite (5.9) with k = k(A) we have instead of c a constant < - 1 as determined 
in the third part of Lemma 5. 19, so that 

dnllf\q)<- Z ^pW) 

m@f.F„(V )3q}{p'£V :F„(p')t<]} 

however on the right the values r g 1 (1 + A) l ^ s+X) and there is at least one point Fgp' from 
V(®k(A))\Vo for which d m {p,F e p') = d kw (p,q)-l, so that ^(Fgp') > C exp(-cd kW+c (p,x)) 
as seen in Corollary 5.21. Since all terms on the right have the same sign, this term gives 
an upper bound. □ 

6. Estimates of the resolvent kernel on the positive real axis 

In this section of the paper we use the estimates of Section 5 and the series expres- 
sion (3.3) for the resolvent G w (x,y) of the Laplacian to obtain estimates in the case 
A € (0, oo). The main result is as follows. 

Theorem 6.1. There are constants k\ and k 2 depending only on the fractal, harmonic 
structure and measure, and such that if A > 0, 

(6.1) (1 + i)- 1/(s+1) exp(-^ l£ / it(/1) (x,3;)) < G (A \x,y) < (1 + A)" 1/(s+1) exp(-K 2 d m (x,yj) 

except that the lower bound is not valid if x or y is in a cell Fg(X), 6 6 ®k(X) such that Fg( Vo) 
contains a point ofVo \ {p}. For these latter cells, the appropriate estimate is instead one 
on the normal derivative. Specifically, if p e Vo, y is not in one of the above cells Fg(X) 
and the normal derivative d' n is taken with respect to the first variable, then 

(6.2) exp(-/Ci4 W )0?,y)) < -d' n G w (p,y) < exp(-K 2 d k(A) (p,y)). 

A symmetrical result holds for the normal derivative —d'^ with respect to the second vari- 
able. There are also bounds appropriate to the case where both points are near the bound- 
ary. If p and q are points in Vq then 
(6.3) 

(1 +A) l/{s+1) exp(- Kl d kw (p,qj) < d';d' n G (A \p,q) < (1 +A) l/{S+1) exp(-K 2 d kw (p,qj). 

It is easy to see these estimates are equivalent to the following global bounds. Let R(Vq, x) 
denote the resistance distance from x to Vq. Then 

exp(- Kl d kw (p,y)) 5 A-V^Rix, V )- l +A l '^)(R(y, Voy'+A^^G^y) < exp(-K 2 d m (p 

With this in hand will be relatively easy to obtain a similar result for the Neumann re- 
solvent Gw(x,y), which has vanishing normal derivatives rather than zero values at points 
ofV . 

Corollary 6.2. For A > the Neumann resolvent satisfies the upper and lower bounds 
of (6.1) everywhere on X. 
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The proofs of these results occupy the rest of this section. We begin by expanding the 
expression (3.3) for the resolvent of the Laplacian and using (3.4) to obtain 

(6.4) G«\x,y)=^ J] r^f^^ x \F^x)^{F^y), 

weW, p,qeVi\V a 

where we recall that G p } is the inverse of the matrix sL defined in (3.5). In the next few 
lemmas we apply the estimates from Section 5 to the terms in this series, for which purpose 
we require the following assumption. 

Assumption 6.3. For the remainder of this section we require that A > 0. 

Lemma 6.4. Let D w be the diagonal matrix with entries D { pp = (1 + /t) 1/(s+1) (B^)~ and 
E (A) = (pMyi _ (1 + /l )-i/(s+i) B (A)_ Then f or a u A and p we have D W ~ 1; E W = 0> and 

for q + p 

e.xp(-e\cd m (p, q)j < Ep\\ < exp[-czd kW (p,q)J, 

where ci and C2 depend only on the fractal and its harmonic structure and c is the constant 
in Theorem 5.1. 

Proof Comparing (3.5) and (3.2) we have 

However Theorem 5.1 then shows that for p — qwe have BpJ (1 + A) l ^ (s+1 \ so Dpp ^ 1. 
Also from Theorem 5. 1 we have for p + q 

-(1 +A)^ s ^exp(-d ku - itJiA) (p,q)) < rfd^fiF-j'iq)) 

< -(1 +A)^ s+l) exp(-cd k(r]fl]A) (p,q)), 

and the observation that d k (j- jtljA )(p,q) — d k ( A ){p,q) lets us choose appropriate constants C[ 
and C2- □ 

Lemma 6.5. For any word w e W* and all A > 0, 

sxp^-cicd k{rwflwA) (p,q)) < ((r w /i w ) _1 + A)~r w G ( pf uA) < exp(-c 2 d k{rwflu , A) (p,q)). 
where C\, C2 and c are as in Lemma 6.4. 

Proof. Recall that d k ( A ){p, q) > e k{X) = (1 + A) llls+l) , so that exp(-d k(A )(p, <?)) can be made 
arbitrarily small by taking A large enough. Lemma 6.4 then implies is close to diago- 
nal, so we may find G w = 1 via the Neumann series. Specifically, 

oo 

(1 + A) 1/(S+1) G W = (/ - D w E w Y l D w = D« + Yp^E^fD^. 

provided A > C\ where C\ is chosen large enough that the series converges. Observe that 
this C\ depends only on the structure of the fractal, because D (A> contains only values ^ 1 
and the values in E w satisfy the estimate in Lemma 6.4. Notice also that all values in D (A) 
and E (A> are positive, hence the same is true of all terms in the series. Making the obvious 
upper and lower estimates of the sum of the series we conclude that for all p and q 

exp(-c lC d k(A) (p,q)) < (1 + /l) 1/(s+1) G^ < exp(-c 2 d k(A) (p,q)). 
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Substituting (r w p w A) in place of A and using r w p w = r^ +1 gives 

exp(-cicd k (r wflw A)(p,q)) < ((r w /i w ) _1 + A)~ r w Gf c f"' A) < exp(-c 2 d k{r ^ A) (p,qi), 

which proves the estimate for A > C\(r w p w ) . However if < A < Ci(r K ,/j„,) _1 then it 
is immediate that ((r„,// vt ) _1 + A)~ r H , is bounded above and below, and we see Gjjf is 
bounded above and below by continuity of its dependence on r w p w A, At the same time 
dk(r„ji w A)(p, q) is bounded because it is the distance on a cellular partition of bounded scale. 
All constants depend only on the fractal and harmonic structure, so the result follows. □ 

Lemma 6.6. There are positive constants C3 and C4 depending only on the harmonic struc- 
ture, such that if c is the constant from Theorem 5. 1 and p e V\\ Vq, then for x in the support 
°f (F~^ x), which is the union of all cells F w j(X) that contain F w (p), 

(6.5) exp(-C3cd k(A) (F K p, x)) < t/^"' 1 ' '(F" 1 x) < e,xp(-c A d k{A) (F w p, xj), 

except that the lower bound does not hold on the cells F w j o Fg(X)for those 6 € ® k (r w -n w ,X) 
such that Fg(X) intersects Vo \ {p}. For these exceptional cells the correct estimate is that 
if x — F w (q)for some q € Vo then 

exp(-c 3 cd m (F w p,F w qj) < -((r^) -1 + Ay lKS+x) d n (^ A) o F~})(x) 

(6.6) ' 

< exp^-c 4 d kw (F w p, F w q)J. 

Proof. From (3.2) we see that \jrp w>1 {F^ x) is a piecewise /l-eigenfunction on F W (X) with 
value 1 at F w (p) and zero at the other points of F W (V\). It is non-zero precisely on the 
the cells F w j(X) such that p e Fj(Vo). On each such cell it is equal to jjWS wMw * /r-j 

According to Theorem 5. 1 we have 
(6.7) 

exp(-cd KrwJtlwjA) (F- l p,F-)(x)j) * ~ ex P(-*(^^)(^7 1 /''^)^)))' 

except that the lower bound does not hold on the cells excluded in the statement of the 
lemma. 

From the partition 0t(r, W r w /< w /i) of X, form the partition r w ® k ( rjf2jruPt[ . A ) of F W {X) and ob- 
serve that the scale is comparable to that of ® k (A.) restricted to F W (X). It follows that there 
are positive C3 and q depending on the harmonic structure so that for y,z £ X, 

(6.8) c 3 d k(A) {F Kj y, F wj z) > d k{ru , jtlujA) (y, z) > c A d k(A) {F wj y, F wj z) 

provided we round appropriately. Substituting (6.8) into (6.7) and eliminating the round- 
ing by taking a suitably large multiple of the exponential gives the desired estimate for 
xj/p^ iF^x) because it is a finite sum of such terms. 

For the normal derivative we note that if x — F w q for q e Vo then x = F WJ q' for some 
q' £ Vq\ {/?}. Using Theorem 5.1 at q' we obtain 

exp(-cd k{rwJ , lrjA) (F^p,q')) < -(1 + r w ^ wj Ay ll ^d n il^\q') 

(6.9) ' 

< exp(-d KrwjflwjA) (F j l p, q')J. 

Recalling that 
we see that 

-((r^jT 1 + ^)- 1/(S+1) 5„(^ ) o F-))(x) 
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also satisfies the estimate in (6.9). As d n \]/p wX ' o F" 1 is a finite sum of such terms, the 
bound (6.6) may be obtained by substituting (6.8) into this estimate. □ 

We divide the proof of Theorem 6. 1 into two parts: the proof of the upper bounds and 
the proof of the lower bounds. 

Proof of Theorem 6.1: Upper bounds. Fix x,y in X. If x + y let w e W* be the longest 
word such that x,y e F W (X), and otherwise let w be an infinite word such that F W (X) = 
M = \y) ( m this latter case the expression (6.10) below may need an additional sum over 
the possible choices of w, but we suppress this because it does not otherwise affect the 
working). Then the series (6.4) terminates at scale \w\ (which is +oo if x = y) and may be 
written 

M 

(6.10) G K, (x,y) = 2_ l 2j r Mi G pq ( F [w), x ^g *■ M,^ - 

i=0 p,<?eVi\V 

It will be convenient to divide the sum into three pieces according to the size of r w p w A. 
Define 



(6.11) i = 



(min{/ : r [w]i p [w]i A < 1} if r H ,p w A < 1 
1 +oo if r K p H ,A > 1 , 



so that if io < oo then r[ w ] ( /I The three pieces of the sum are as follows, where 

we note that it is possible for a piece to be empty. 



brl 



h = 



min{io— 

A L,.,-T7 \ T/ 



o p^ € yi\vb 

min{io-l,|w-i} 



/ 3 = 



A „,-17 \ T7 



/>eVi\Vo 

For both l\ and /2 we use the trivial estimate 



for the i// factors. Then for /] we have r^yp^A < 1, from which the values |G 
1 , and so 



,(r lw] .n lw]t X)i 
PI 



M 



(6. 12) < /, < J] r Mi < r MiQ < (1 + A)" 1 '^ if r w p w A < 1 

i'o 

/i = if r w p w A > 1 

where we used the geometric decay of the r^ w y and the definition of i'o. Note that we have 
the upper bound (1 +Xf x ^ s+V> rather than just A~ lli - s + l) because when A < 1 we have iq = 
and r[ w j. = 1 corresponding to the empty word. 
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For I2 we instead estimate the factors r[ w ]. q^ 1 *™'* using Lemma 6.5 to obtain 

min|! -l,M} 



(6.13) 



2 2 r[w] ' G " 

ptqeVAVa 
min(! -l,|w|) 

^ Yj (Om./Vl)" 1 + i)" 1/(S+1) exp(-C2d*( WM ^)(p. ?)) 



f(l + A)" 1/(S+1) ifr w /i w A<l 
[(1 + exp(-c 2 d*( rw ^)(/>, qj) if r wj u w /l > 1 

because d^p, q) is at least exponentially increasing in k, so the sum is dominated by a con- 
stant multiple of the largest term, which is the one with the maximal value of ;. Moreover 
the constant multiple depends only on the r and p values and the way path lengths on &k 
grow with k, all of which are properties only of the fractal, harmonic structure and mea- 
sure. In the case that the maximal value of i is ;'o we have r w p w A =* 1 so d^ rwflw X)(p, q) < 1 
and the exponential term is trivial. The factors (1 + AY x ^ s+l) can be used rather than the 
more obvious choice A~ x ^ s+l) because if A < 1 the sum is empty. 

For the ^3 term we must take a different approach, because in this case i < io implies we 
can apply Lemma 6.5 to obtain 

and therefore the only estimate we have is that from Lemma 6.6, which gives 

min|/ -l,M) 

/ 3 < (1 + Ar^ 2 2 #*^(P[iD ^ Wi Vf>) 

i=0 peVAVo 
min{i -l,M) 

(6.14) <(l+^)- 1/(5+1) J] J] exp(-C4^ w (F [H , ] ,/ 9 ,x)-c 4 4(, 1 )(^ M ,P,3')) 

i=0 peVAVo 

where we again have replaced the obvious A ^ s+ ^ factor with (1 + A)~ l t (s+1) because the 
sum is empty if A < I. 

The above indicates that we need a lower bound on d^)(F[ K ].p, x) + dk(x)(F[w],p,y) for 
the points F^p, p e V\ \ Vo and a mechanism for counting how many of them there are. 
Let Li be the set of such points that are in the boundary of the /-cell containing x but not 
in the boundary of the (/ + l)-cell containing x, then the number of points in L; is smaller 
than the number of points in V\ \ Vq, which is a constant depending only on the structure 
of the fractal. Moreover any point z e L; is separated from x by an (/ + l)-cell. By the same 
reasoning as in Lemma 6.6, the d^ X ) diameter of any such (/ + 1) cell F[ W ] U -(X) is bounded 
below by a constant multiple of the dk(r lw} p lw] x> diameter of X, rounded down to the nearest 
integer. Writing the diameter of X with respect to dk as diam,t(X) we obtain: 

Yj vw{-cd K A)(F m p, x)) < exp(-c dmm kirMiMMiA) (X)). 

We use this estimate only for I < \w\. Any points F [ W ].p which occur in the sum but are 
not in any of the L/, / < \w\ must be in F W (V\), so the number of these is bounded by a 
constant depending only on the structure of the fractal. For these points we use the triangle 
inequality to estimate the corresponding terms of (6.14). 

exp(-C4d kw (F lw]i p,x) - c 4 d kW (F [w]i p,y)) < exp(-c 4 ^ (/1 )(x,y)). 
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(6.15) 



Combining these estimates and substituting into (6.14) we have shown 

min{i ,M} 

7? < (1 + Ay 1/(s+l} (exp(-c 4 d k( A)(x,y)) + ^ exp(-cdiam Wr[u , |(jU[u , ]( . 1) (X))j 

1=0 

(1 + ^-l/cs+i) ifr w /z w A<l 
(1 + A)-V( s+1 ) exp(-c 4 <4 W) (*,y)) if r w fi w A > 1 

where we used the same reasoning about the exponential decay as was used for the h 
term, along with the fact that d k (A)(x,y) < diam^.^^X) because x and y have resistance 
separation < r w . 

Comparing (6.12), (6.13) and (6.15) we see that we always have a bound by (1 + 
^-i/os+i) jf fynfi^A < 1, and in this case the resistance separation of x and y is at most 
a constant multiple of r w = (r H .. / u M ,) 1/(s+1) < A~ 1/(s+l \ so that d k(/X) (x, y) < 1. Setting 
K2 = min{c2, C4} we therefore obtain 

CP>(x.y) < (l + A)-^ s ^ < (l + Ar^e W (-K 2 d kw (x,y)). 

For the case r w fi w A > 1 we have I\ = 0, and the estimate for I2 is dominated by a multiple 
of that for 73, so again 

G w (x,y) < (1 + Ar 1/( - s+1 ^xp(-K 2 d k(A) (x,y)) 

and we have established the upper bound of (6. 1) stated in the theorem. 

We now turn to the upper bounds for normal derivatives. Recall that d' n denotes the 
normal derivative with respect to the first variable. If x e Vq and y + x then the series 
for G w (x,y) is finite and we can compute the normal derivative term by term. The upper 
bound on -d' n G^ A \x, y) may then be obtained in almost the same way as the upper bound on 
G w (x,y). The reasoning for both h and 73 is unchanged except that instead of bounding 
r Mi Gp™' /Vli/i) by we use the full upper bound ((r^j./^].)- 1 + Ay 1/(S+1) from 

Lemma 6.5. The upper bound for -d'^p" 1 '^'^ ° FT^ (x) from (6.6) cancels this factor 
and leaves precisely the exponential decay terms seen in (6.13) and (6.14), so the only 
change to these estimates is that the / t _1 /< s+1 > factor is no longer present. The 7i term 
requires slightly more changes. Again we use Lemma 6.5 and (6.6) to see that 

(6.16) r [wL G^ IiMMiA) (-5>^ MiMM ^ ) o F-^( X )) < cxp(-c 4 d m (F [w]i p, x)) 

is bounded, but then use it to conclude that 

M 

Kh<J^ J] exp(-C44 W )(f Mi P^))^ r[ " I, " Mi ' l) (i 7 M 1 j) 

i'o p,qeVi\V a 

< exp{-c A d m {F[w],P,y) ~ c 4 d k( A){F [w]i p,x)), 

where the last step is from (6.5). This sum can be bounded by the same argument as was 
used for ^3 in passing from (6.14) to (6.15), so we may sum the terms from /1, I2 and 
73 to obtain the desired bound -d' n G w (x, y) < exp(-K2<4(i)(*, y) ) provided y + x. This 
verifies the upper bound in (6.2). An easy argument then shows that for points z within 
the ® k (A) cell containing x, the peak size of G w (z, •) is comparable to R(x,z) rather than 
(1 + Ay l t (s+l) ; we will later need the immediate consequence 

(6.17) \\G w (z,-)\\v <(l+A)- 1/ ( s+l) R(x,z)asz^ xeV . 
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To complete our upper bounds we must deal with what happens when x and y are both 
in Vq. If this is true and they are distinct, then the sum (6.10) reduces to the finite number 
of terms with i = 0. We can break each term into the product of a factor like that in (6.16) 
and a factor for the normal derivative in the second variable, which is of the form 

-d:;(^ mA) o F^)(y) < ((r^MM,)- 1 + A) 1/{S+V) exp(-c 4 d k(A) (F lw]i p,y)) 

by (6.6), but with i = 0. However, the points F[ W ] p = p are then in V\ \ Vq, so both 
dk{A){F[w]iP,y) and dk(,\){F[ W ] p, x) are comparable to d^ix, y), and as ; = the overall 
bound reduces to 

d' n d';,G w {x,y) < (1 + A) 1/(s+n exp(-c 4 d kw (x,y)). 

Finally, we look at d'^d' n G w (x, x) for x e Vq. Notice that in the above discussion of the 
normal derivative d' n G \x,y) we established that (summing over those w corresponding to 
x if necessary) 
(6.18) 

A'r^W, \ V V r> ( W»li ,t ) / of , (l-lirt"'!," 1 ) 17-1 i \\ ,( r Mi/Vi i ' i )/T7-l s 

- d n G w (x,y) = 2_, £j r Wt G Pi \-°r$p ° F M i Wj^ r ? ( p iw] t y) 

i p,qeV,\V 

in which the sum is finite for any y + x. However replacing ^^"^"^ (F^\ y) with 

— d^iftq™ 1 ^' (.FT^i y) throughout the series does not give a series for d'^d' n G (X \x,y). The 
reason is that every finite truncation of the series is zero at x, whereas -d' n G (A \x, x) is the 
continuous extension of -d' n G w (x, •) from X \ \x] to x, which exists because the series is 
a sum of continuous positive functions in which all partial sums are < e.xp[—cndk(A){x,y)) 
for y + x. We therefore compute d" } d' n G (?C) {x, x) using 

&'d'p x \x,x) 



(6.19) 
(6.20) 



where hj is harmonic on F[ W ].(X) and equal to 1 at x and zero at the other points of F[ W ].(Vq). 
We used that 2vef m (v ) d„h*(y) = because h* is harmonic, and that d' n G w (x, x) = 
\m\ Zr ^ x d' n G iA \x,z) by definition. Now rewrite these terms using the Gauss-Green formula 
in the first variable. Since (A - A')G w (s,y) = 6 v (s) we have 

& n G w {x,y)= f Us)A'G w (s,y)dp.(s) = A f U^)G w (s,y)dp( s ) - ( x (y) 
Jx Jx 



= limE(-d'„G {A \x,-),h*(-)) 

i—*oo 

= lim V d n h;:(y)(-d' n G w (x,y)) 

I— iOO ' * 

= lim V d„h*(yMG w (x, x) - d' n G w (x, y)) 

I— >oo < * 

= lim lim V d n h*(y)(d'„G w (x, z) - d' n G w {x, y)) 



yeF M .(V )\{x} 
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from which (6.20) may be written 

d:;d:,G w { X ,x) 

= lim V d n h*(y)(t x (y)-]imt x (zj) 

y£F lw]i (V )\{x) 

-Alim f Us) X d„hJ(y)G (A \s,y)dfi(s) 

ux yeF M .(V )\f;t} 

+ A\im y d n h*(y)]im f Us)G w (s,z)d^s). 

yeF Mi (V )\[x\ 

(6.21) = d„Ux) - A lim f V d n hf(y)G w (s,y)dKs) 

l— >oo l v X i 

yeF m (V )\{x) 

In the last step of this computation we used that lim-^ x Uz) = U x )> so the first term just 
produces d„U x ) by the same argument that gave (6. 19). We also found that the third term 
was zero, because of (6.17) and the boundedness of h*. 

To further simplify (6.21) we again apply the Gauss-Green formula, this time in the 
second variable, to obtain 

J] d n h*{y)G {A \s,y) = d';G (A Xs,x)- f h*(t)A"G w (sj)dKt) 

= d" l G (A \s,x) + h>;(s)-A f h x i (t)G {A \s,i)dn(i) 

Jf M jQ0 

Two pieces of this make no contribution in the limit as i — » oo. One is the integral 

A I ^(i)/i; c (i) c/yu(i) <K F M(X)) -» as i -> oo. 
Jx 

For the other we use the upper bound of (6.1) to compute f G^is, t)dfi(s) < AT 1 and 
therefore 

A 2 f Us) f h*(t)G w (s,t)dn(t)dKs)lZA f ^)#)<^(F W( x)) 

JX JF Mj {X) JF Mj (X) 

which also goes to as i — > oo. We used that hj is bounded by 1. On the remaining piece 
we use the previously established upper bound in (6.2) to see that that — d'^G^'(x,y) has 
faster than exponential decay with characteristic length k(A), so that its integral against a 
bounded function is estimated by the integral over a cell of diameter k{A) at x, which has 
measure comparable to A~ s ^ s+l \ Thus 

(6.22) A f Us)(-d',;G w (s,x))dfi(s) < A l ' iS+l \ 

Jx 

Substituting all of this into (6.21) we have at last 



(6.23) dy„G«\x, x) = d„Ux) + A- f Us)(-d';G w (s, x)) dn(s) < (1 + A) 1 

Jx 



-,1/(5 + 1) 

Ix 

This completes all the upper bounds in the theorem. 



Proof of Theorem 6.1: Lower bounds. Observe that all terms in the sum (6.10) are posi- 
tive, so to obtain the lower bound of (6.1) it suffices that we have the desired bound on 
a single term. Our upper estimates suggest we consider the dominant term for both the 
diagonal and off-diagonal series. To do so we will have to split into two different ranges 
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of A. This splitting depends on the lower estimate from (6.5) of Lemma 6.6, which is only 
valid off some exceptional cells. 

The exceptional cells in the sense of Lemma 6.6 are of the form F[ W ].j o Fg(X) for those 
6 e ®k(r Mjjt j Mjl A) such that F e (X) intersects Vo \ {p}- Observe that if r^jfi^jA is too 
small then these cells could cover F[ W ] t (X), so to use this estimate at all we must assume 
that r[ W ]./j,[ W ].A > c, where c is chosen large enough this does not occur, and is a constant 
depending only on the fractal and harmonic structure. 

Fix two points x and y and recall w is the longest word so F W {X) contains both x and y. 
Our two cases are that where r w p w A > c and where r w p w A < c, the latter including the case 
x = y. In the former we prove the lower bound (6. 1) including the exponential term, but in 
the latter the exponential term is bounded below and plays no role in the estimate. 

First consider the case r w fX w A > c. Then the lower estimate from (6.5) of Lemma 6.6 
is valid off the exceptional cells of F W (X), and these cells have size comparable to cells of 
the partition ® k ( A ). Fix such a size and call the subcells of F W (X) that both have this size 
and contain a point of F w (Vq) the boundary cells of F K (X). We further arrange this size so 
that the exceptional cells for ij/p^ ° F~^ are contained in those boundary cells that do not 
contain F w p. 

From Lemma 6.5 and Lemma 6.6 we know that for any i such that r[ w ] t M[w]i^ > c and x 
and y not in the exceptional subcells for tfr^' j Mi o 

Z ((Vl^m, + AY U(S+1) exp[-Cid k (r M .n lw]t A)(x,y) - c 3 cd kw (F [w]i p, x) - c 3 cd k(A) (F [w]i q,yi) 
(6.24) 

> A~ ll(S+l> exp(-cic 3 cc/i ( , 1) (F M ./?, />],<?) _ Cicd k ( X) {F[ w] .p, x) - c 3 cd k{A) (F [w]i q,y)j 

where the dk(r lw] .n lwi .A.)(x,y) term from Lemma 6.5 was converted to c^cd^iF^p, F^q) 
using (6.8). To get our lower bound (6. 1) on G iA \x, y) for A > 1 it suffices that given x and 
y, neither of which is in an exceptional cell containing a point of Vo, we can find such i, p 
and q with x and y not in exceptional subcells of F[ W ^ (X) and 

(6.25) d m (F [w] .p, F [w] .q) + d k(A) (F [w]i p, x) + d k(A) (F [w]i q, y) < d m {x, y) + k' 

for some constant k' depending only on the fractal and harmonic structure, as the corre- 
sponding term from (6.24) will be a lower bound of the correct form for the non-negative 
series (6.10) for G w (x,y) once we set k\ = c\Ct,c. 

One possibility is that neither x nor y is in a boundary cell of F W (X). Our choice of 
w ensures that x and y are in separate level-one subcells of F W (X) (i.e. cells of the form 
F w j(X)), so the d k ( A ) geodesic between them must pass through a vertex F w p for some 
p e V\ \ Vo- We conclude that d k ( A )(F w p,x) + d kiA) {F w p,y) = d k ( A )(x,y), which is the 
desired bound (6.25) in the case q — p and i = \w\. 

The alternative is that either x or y or both is in a boundary subcell of F n (X). There is no 
loss of generality in assuming x is in a boundary subcell, and the corresponding boundary 
point is F w p' for some p' e Vo- We will assume that F w p' £ Vo, in which case there is a 
largest ; satisfying F K p' — F[ W ].p for some p e V\ \ Vq. If y is also in a boundary subcell 
then by taking x to be the point with the larger value of i of this type, we can assume that 
the boundary point corresponding to y is of the form F[ W ] t q for some q e V* \ Vo- 

Since the boundary subcells have size comparable to cells from ® k ( A ) they have bounded 
diameter in the d k ( A ) distance. Therefore d k ( A ){F[ W ].p, x) < ko and the triangle inequality 
ensures d k ( A ){F[ W ].p,y) < d k ( A )(x,y) + ko. Examining the term of (6.24) with this i and p 
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and taking q — p we find that x is too close to F^p to be in an exceptional subcell for 
^.wjiMM; o If y is not in an exceptional subcell for this function then 

dkcX)(F [w] .p, F M .p) + d k{A) (F [w] .p, x) + d kw (F[«],p,y) < k + d k( A)(x,y) + k 

which verifies the desired bound (6.25). However for y to be in an exceptional subcell 
for i/fp^ wV o F^. it must be in a boundary subcell for F[ W ] ( (X) and thus for F n ,(X). We 
conclude that i + 1 = \w\ and therefore y is within d^x) distance ko of F[ W ].q for some 
q e V\ \ Vq. Looking at the term from (6.24) with this i, p and q we use the triangle 
inequality d kW (F M .p, F [w]j q) < d k(A) (x,y) + 2ko to obtain 

(6.26) d k (A)(F lw]l p, F M .q) + d k (A)(F [w]i p, x) + d k (A)(F [w]i q, y) < d k (A)(x, y) + 4ko 

which verifies (6.25). Then x is not in an exceptional subcell for ^* 1 ^ W(/1 Fj~^. and y is 

not in an exceptional subcell for ^ wli/iWi Fj"^ , so the lower bound follows from (6.26). 

Thus far we have verified the lower bound of (6. 1) in the case that r w p w A > c and under 
the assumption that if either x or y is in a boundary subcell of F U .(X), then this subcell does 
not contain a point of Vq. However in the case r w p w A < c the lower bound amounts to 
saying that G (,t) (x, y) is bounded below if neither x nor y is within a 0*^) cell containing a 
point of Vo, which is vacuous when A is so small that these cells cover X, and evident from 
continuity of G w {x,y) on the bounded interval in A otherwise. 

Next we examine what happens on boundary subcells which contain points from Vo- In 
this case we seek lower bounds for the negatives of partial derivatives as in (6.2) and (6.3). 
Again we assume r K p w A > c so as to use the lower bounds that go into (6.24). 

Recall that if x e Vo and y + x we can write -d' n G^ } (x,y) as the series in (6.18). For 
this sum we can use termwise estimates like (6.24), because replacing the lower bound on 
^MiWH-i (jCj-i x ) in (6.24) by the lower bound on its normal derivative from (6.6) produces 
the same result except without the factor A~ l ^ s+l \ As in the reasoning leading to (6.25), 
given x e Vo and y not in an exceptional cell containing a point of Vo, our lower bound 
on —d' n G w (x,y) will follow if we can find i, p, and q with y not in an exceptional subcell 
of F[ W y(X) and such that (6.25) holds. Moreover the situation is simpler than it was there, 
because x will always be in an exceptional subcell of F^iX) and we can take p so F[ W ].p = 
x. For the largest number j such that y 6 F[ W ] (X), if y is in an exceptional subcell of 
F[ W ]j(X) then it is not in an exceptional subcell of F[ W ], , (X), so we may take i to be either 
j or j - 1 as appropriate. In either case it is apparent that (6.25) is valid, establishing 
the lower estimate of (6.2) when r w p w A > c. To deal with the case r w p w A < c we fix 
y + x and observe that G (x \y,y) converges to a positive constant as A — > and that 
(-d'„G^ l \x, y))(G ( ' !) (y, y)) is positive and decreasing in A by an argument similar to that in 
Corollary 5.8. An easy computation using the series (6.18) with A = shows -d'„G <0 \x, y) 
is positive, and we conclude that -d' t1 G w (x,y) is bounded below by a constant depending 
only on the fractal and harmonic structure. Together with the above, this shows that the 
lower estimate of (6.2) is valid for all positive A. 

We can now use this lower bound to improve the < in (6.22) to ^, because it guarantees 
that —d'^G^ix^) achieves size comparable to a constant depending only on the fractal and 
harmonic structure. The lower estimate on the integral then follows by an argument like 
that in Lemma 5.9. However this and (6.23) immediately imply the lower bound in (6.3), 
so all lower bounds in the theorem have been proved. □ 
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Proof of Corollary 6.2. According to Theorem 4.2 of [6], the Neumann resolvent is ob- 
tained from the Dirichlet resolvent via 

G%\x,y) = G«\x,y) + £ C^f(xWf(y) 

where Cp) are the entries of the matrix inverse to that with entries d n r/^\p). From The- 
orem 5.1 the values cL are all positive and the diagonal values C p p are comparable to 
(1 + A.y ll( - s+l \ From the same theorem it is evident that J] ( p A> (x)T] ( ^\y) is bounded above by 
exp(-K2<4(,i)(X;y)), so the terms added in passing from G w (x,y) to G^\x,y) also satisfy 
the upper bound (6. 1), and this estimate is therefore valid for G\t (x, y). As we are adding 
positive terms, the lower bound for G w (x, y) in (6. 1 ) applies also to G^ (x, y). To complete 
the proof we need only check that this bound is still valid for x and y within a cell of scale 
k(A) at a boundary point. At such points the exponential term is trivial, so the estimate is 
just that G£ (x,y) > (1 + X)~ l ^ s+r> . This is true because rfp is comparable to 1 at both x 
and y from the estimates in Theorem 5.1, andCp'^ - (1 +A)~ 1 ^ S+ ^ as previously noted. □ 

7. Resolvents on Blowups 

Definition 7.1. For an infinite word o~ = 0"i0"2 • • ■ let 

G-n = (F-% ]n (X) = F~\ o ■ ■ ■ o F~ l n (X). 

We call Q_„ a finite blowup of X. Observe that the sequence {Q.-,,} is increasing because 
Fj(X) c X for any j. The blowup Q. of X is the union 

a = (J q_„. 

n>0 

In general there are uncountably many non-isometric blowups of a fractal. If the word 
ends with the infinite repetition of a single letter, then the blowup will have non-empty 
boundary, but otherwise there is no boundary. One reason to be interested in blowups is 
illustrated by the familiar case of the unit interval [0, 1], which is a pcfss invariant set of 
the iterated function system Fq(x) - |, F\(x) = ^j- on R. In this case any blowup corre- 
sponding to a word that terminates with 000 ■ ■ ■ is a ray to +oo, while if the word terminates 
with 1 1 1 ■ ■ ■ it is a ray to -oo; for all other words the blowup is all of R. Heuristically we 
think of the relationship between a fractal and its blowup as somewhat akin to the relation- 
ship between [0, 1] and R, inasmuch as they are compact and non-compact examples of 
sets with similar local structure. Some interesting results about blowups of fractals may be 
found in [17, 20]. 

Fix an infinite word cr. It is straightforward to define a Laplacian on the blowup Q 
that is consistent with that on X. Recall that for a finite word w, we have A(m o F w ) = 
r w fi w (Au)o F w . Now if u is defined on £2_„ then u °(F~ l )[o-] n is defined onX, so to maintain 
this composition rule for the Laplacian we must set 

Aw = r Mn/U[(r]n (A( M o (F" 1 )^)) o ((F" 1 )^)" 1 

(7.1) = JVlWo-],, F<T„ ° Fo-^ O ■ ■ ■ O Fo-j . 

It is easy to check that this ensures the definition of A on fi_„ is consistent with that on any 
other Q_ m , so suffices as a definition of A on O which is consistent with that on X = £2o- 
We may similarly define E and ponQ. (note that then is infinite but cr-finite) so as to be 
consistent with their definitions on X. 
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Now the result (3.3) from [6] may be transferred to each £!_„. For convenience of 
notation, for the remainder of this section the function G w (x,y) of (3.3) will be renamed 
Gq (x,y) to emphasize that it corresponds to the zero blowup X = £2o- 

Lemma 7.2. Suppose z e C is such that there is no finite word wfor which (?"[o-]„^[o-]„ )~ 1 r w p w z 
is a Dirichlet eigenvalue of A on X. Then defining 

(7.2) G { *(x,y) = r^G^'^iF^ ° • • • o F^x, F an o • • • o F aJ ) 
the function 

u(x) = f G ( = ) „(x,y)f(y)dp(y) 
Jn_„ 

solves (zl — A)u — f on Q „ with Dirichlet conditions at 5((F _1 )[o-] n (X)) = (F )[(,-]_ (Vo)- 
Proof. By direct computation setting xf = F^ o ■ ■ ■ o F tTl x and similarly for y, 

(zl - A x )u{x) = {zl - A x )r^ f G^'^V^ ° • • • o F^x, F an o ■ ■ ■ o F ai y)f(y) dp(y) 

Jn_„ 

= (2/ - AJr^-J! L f G^ h Vcr, o ■ ■ • o F^x,/)/^ 1 )^),/) <W) 
- (r^j^-A,) f G^^^x'.yO/CCF-^M^O^Cy') 

= /((F _1 )[tr]„^') = fix), 

because (r^\ p.^\ z) satisfies the hypotheses of Proposition 3.2. The Dirichlet boundary 
conditions are immediate from the original result. □ 

Theorem 7.3. For any A > the sequence G_„(x,y) is uniformly convergent on £1 x CI. 
The limit G^\x,y) is such that if f e L'(Q, dp) then 

«(*)= f G^{x,y)f(y)ditSy) 
Jo. 

satisfies (AI — A)u — f on Q. IfQ, has boundary points then u — Oat these points. Moreover 
there are positive constants kt, and K4 depending only on the fractal and harmonic structure 
and such that 

(7.3) A- 1/(S+r, exp(-K 3 d kw (x,yj) < G ( *\x,y) < AT ms+1) exp(-K 4 d m (x,yj). 

unless x or y is in a cell of & k (A) that contains a boundary point. If there is a boundary 
point p and y is not in such a cell, then we have 

(7.4) exp(-K 3 d m (p, yfj < -d'„G^\p, y) < exp(-K A d m (p, y)) 
and if p,q are boundary points then 

(7.5) exp(-K 3 d kw (p, q)) < W n G { *\p, q) < A ll{s+X) exp(-K 4 d kw (p, q)). 

Proof. We need only establish the convergence and the bounds; the other assertions follow 
immediately from Lemma 7.2, which is applicable because A > and the spectrum of 
A contains only negative values. Observe that instead of x,y in a general compact set it 
suffices to consider x, y e X = £lo, because the proof is the same for any Q_„ and any 
compact set is contained in some sufficiently large Q_„. 
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For the moment we assume the convergence and examine the bounds. These are readily 
deduced from the estimates in Theorem 6. 1 applied to (7.2). Direct substitution into (6.1) 
shows that both the upper and lower bound for G_l(x,y) are of the form 

'm/ 1 + r M»^M»' 1 )" 1/(S+1) e M~ Kd «ti h A^A Fcr " ° " ' ° F(T ^ F,T " ° " ' ° F<T ^) 
however 

-l /i , -l -l n-W+l) r , ,x-i/(s+D i-i/(S+i) 

and 

(7-6) d k(r-^ hi A)( F <r„ ° ■ ■ ■ ° F^x, Fa- n ° ■ ■ ■ ° F^y) * d m {x,y) 

as described in the proof of (6.8). This demonstrates (7.3). 

Now using the fact that the normal derivative of a function composed with Fa-,, ° ■ ■ -d-Tv, 
is the normal derivative of the function times the factor we see from (6.2) that the upper 
and lower bound for — d' n G_„(p,y) are both of the form 

exp(-/a4 ( ,.^^ |n/1) (/v„ o ■ ■ ■ o F ai x, F lTn o ■ ■ ■ o F^y)) 

so (7.4) follows using (7.6). Using the scaling of the normal derivative twice we have 
from (6.3) that the upper and lower bounds for d'^d'fi^ip, p) are both of the form 

n^S 1 + r M„/ i M^) 1/(S+1)eX p(~^'t('-fJ! ] „< ] ^)( F ^ ° " ' ° F <riX> F <r„ ° ■ ■ ■ ° Fo-ji) 
so using 

r lo-],A l + r [a-] n f X [cr]„ A ) ~ Vlo-lMa-l, + *) ~* A 

as n — » oo and (7.6) again we get (7.5). Note that these latter estimates are only of interest 
if there is at least one fixed point p that is in all but finitely many Q.-„, from which it follows 
that p is a boundary point of Q. 

In order to establish the convergence it helps to reorganize the notation so that it is clear 
what we are summing. From (7.2) and (3.3), 

G^(x,y) = rj^G^^Vo-, ° ■ ■ • ° F lTl x, ° ■ ■ ■ ° F^y) 

= rj^ J] r/*^;; 1 o Fo-, o ■ ■ • o F cri x,F-' o F^ o • •• o F^y) 

weW, 

however if w has length m then F~ l o F (J - n o ■ • • o Fq-, = F" 1 o ■ • • o F~J o F ail o ■ • • o Fq-, , 
so produces a non-trivial term in the sum only if w begins with <x„, cr n cr n -\, and so on. 
More precisely, if we think of <x as infinite to the left, cr - ■ ■ ■ cr^crxcri, and write [cr]_„ = 
cr„ ■ ■ ■ cr\ and [cr]_„ „, = cr„ ■ ■ ■ cr„_ m+1 then the words producing a non-trivial term are of 
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the form [cr]_ n>m for < m < n, or [cr]_„w for some w e W*. Consequently 

G W „(x,y) 

r£ u £ Vl- <PW^ ^-^>((F [(r] _ )-i o F M _„*, (F M _ )"> o F M _„y) 

=o 

2 ^-^^-"^-""^-'^-"^((^m^w)" 1 o (F^)" 1 o F [o--]_ n y) 



m=0 

-1 



m=0 weW, 



m=—n weW* 

In particular, for «' > n we have from (3.4) 

-ii-i 

G w nl {x,y)-G w n {x,y) = £ ^fi^^(F w/ ,F w ,j) 



m=—n' 

-n-\ 



m=—n' p.qeViXVo 



From Lemma 6.5, 



and applying Lemma 6.6 we find 

\G^(x,y)-G^(x,y)\<A-^ £ £ ^^V^^f^Vw^) 

m=-n' p,^eV 1 \V 
-n-1 

(7-7) <^ I/(5+1) 2 J] exp(-c 4 4 w (Ff; ]m p,x)- C 4*(i)(^rA^^))- 

m=—n' p,qeV]\Vo 

Recall that we were able to assume x,y e Qo. Now we wish to count how many FjV /? 
and Fri| g are within some distance of Qo. It is not easy to do this using the dk{x> distance, 
but it is easy to make an estimate using the number of cells that have the same scale as Qo, 
by which we mean the cells of the form Fj"^ o F„,(Qo) where \w\ = m. For cells of the 
same scale as Qo, we think of those that intersect Qo as forming an annulus of size 1, those 
intersecting this annulus (but not Qo itself) as forming an annulus size 2, and so forth. It 
is apparent that the number of points in the intersection of {F^ p : p e V\ \ Vo,m > 1} 
with any such annulus is bounded by a constant depending only on the fractal. We can also 
obtain a crude estimate on the growth of the dk(A) distance from Qo to the « th annulus, just 
by noting that the resistance distance to this annulus must grow at least like (max ; rj)~ n , 
and the dk<j)(x,y) distance is bounded below by (1 + A) l ^ s+l) times the resistance distance, 
so also grows geometrically. It follows that (7.7) converges (pointwise) as n' — > oo and 
satisfies a bound 

-c 4 min{^ ( ,i)(F [o : ]m p,x) : m > n, p e V\ \ Vo}j. 
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However we established above that there are only finitely many points FT^ p within any of 
our annuli. Thus for any prescribed distance we can take the finite union of annuli covering 
that distance around Oq, and by taking n large enough we can be sure no FjV p with m > n 

lies in this union. Hence G_„(x,y) converges uniformly to (x,y). □ 



8. PHRAGMEN-LINDELOF TYPE THEOREMS 

In this section we prove some complex analytic estimates related to the Phragmen- 
Lindelof theorem. Later these will be used to prove bounds for the resolvent in a sector of 
C that omits the negative real axis. Phragmen-Lindelof theorems have previously been used 
to obtain off-diagonal decay estimates for heat kernels from Davies-Gaffney estimates [2], 
however it seems that the techniques and results proved there are not applicable in our 
situation. 

Fix an angle a < n and consider the sector 

A a = {z = \z\e ip eC:-a<p<a). 

It will also be convenient to identify A+ = A a n {3z > 0} and A~ = A a n {3z < 0). 

One version of the classical Phragmen-Lindelof theorem is as follows; a proof may be 
found in [14]. 

Theorem 8.1 (Phragmen-Lindelof). Suppose that u(z) is a function analytic in an open 
sector of angular size a and continuous on the closure of the sector. If u is bounded by M 
on the sides of the sector and satisfies \u(z)\ < C exp(c\z\ a ) for some a' < 2 and constants 
C and c, then \u\ < M on the closed sector. 

As a particular consequence we see that 

Corollary 8.2. Suppose g(z) is analytic on A* and continuous and bounded by 1 on the 
closure. If\g(z)\ < exp(— fli|z| a2 ) on the positive real axis for some < a2 < I and a\ > 0, 
then for f5 e [0, a] we have 

|§(|z|eHl < exp — \z\ ai 

\ sin(fl2a') / 

Proof. With the convention that all power functions are defined by cutting the z-plane 
along the negative real axis, let 



ja\ 

v{z) = exp — 
V f 



-z a 



sin(fl2<3') 

which is analytic in A+. Note that the real part of the exponent vanishes on Arg(z) = or, and 
that on the positive real axis it is fli|z| fl2 . Thus |v(z)g(z)| < 1 on <9A+. It is also apparent that 



|v(z)| < exp( , Qx \z\ m \ 



\ sin(fl2<3') 

within A J. Since «2 < 1 and the sector has angle a < n, the Phragmen-Lindelof theorem 
implies |v(z)g(z)| < 1 on the closure of A+. Thus for z = \z\e'P, < f3 < a we have 

i f i t m-1 /-aisin(a 2 (ar -/?)) \ 

\g{z)\ < |v(z)| = exp — |z| fl - 

V sin(a 2 a) / 
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We shall later see that this could be applied to both piecewise eigenfunctions and the 
resolvent in the case of affine nested fractals, because these have decay exp(-aidk ( A)(x,y)) 
for A on the positive real axis and dk(A){x,y) - (1 + A) y ^ s+l 'R(x,y) y by Proposition 4.3 and 
Definition 5.12. In the more general pcfss case it cannot be used because we do not know 
whether dhX)(x,y) grows like a power of A. In order to deal with the more general decay 
exp(-fl 1 c/fc / i)(x,y)) we need a modified version of Corollary 8.2. The modification takes 
into account the fact that dk { A){x,y) behaves like a power of A over exponential scales, but 
that the powers could perhaps be different for distinct scales. In essence, what we must 
do is create functions analytic on a sector and growing according to different powers on 
different exponential scales. This is done using Schwarz-Christoffel functions. The author 
anticipates that the existence of maps might be known, but does not know of a place where 
they are described. Our replacement for Corollary 8.2 in the general case is as follows. 

Theorem 8.3. Fix a < n. Suppose f : [0, oo) — » [0, oo) has the property that there are 
< j8i < /?2 < 1 and C\ > and c 2 > such that 

(8.1) C1 M^^, C2 M*. 

for any sufficiently large M 6 [l,°o) and all j € N U {0}. If g(z) is a function that is 
analytic on A*, continuous and bounded by 1 on the closure of A^, and satisfies \g(z)\ < 
exp(— f(\z\"' a )) on the positive real axis, then there are constants c and c such that for 
y8e [0,a] 

\g(\z\e il3 )\ < exp(-c/(| Z r /a )sinc(l - ^)) 

Proof. From Lemma 8.5 below we know that under these hypotheses for every < e < 
| minjC/?!), (1 — /S2), O62 - fix)} there is a function F(z) analytic in the upper half plane and 
satisfying 

c 3 (e)f(\z\) sin^ - e)(n-/3)) < KF(\z\e il3 ) < /(|z|) sin(08 2 + e)(n-/3)). 

Consider 

( . ( F(z n/a ) \ 



i(H0 2 + e))' 

which is analytic in A*. Observe that z"^ a takes the ray with argument a to the negative real 
axis, so that F(z"^ a ) is imaginary and |v(z)g(z)| < 1 on this ray. Since the real part of F(z K ^ a ) 
is bounded by f(\z\"^ a ) sin(n(fi2 + e)) on the positive real axis, we see that |v(z)g(z)| < 1 on 
dAt. We also have the bound 



|v(z)| < exp . 

V sm(n(J3 2 + e)) / 



on A*. Since f(\z\" /a ) < c 2 \zf 2nla for all sufficiently large |z| and ^ < 2, the Phragmen- 
Lindleof theorem implies that |v(z)g(z)| < 1 on A+. We conclude that 

, , „ , , /-c 3 (e)/(|zr /ff )sin((/? 1 - e)n{\ - 

\g(z)\ < |v(z)| 1 < exp — — 

\ sm{K(J32 + £)) I 

and note that the constants are non-zero. □ 



The construction and properties of the Schwarz-Christoffel function used in the proof 
are contained in the following lemma. 
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Lemma 8.4. Given < a\ < a 2 < 1 and numbers tj such that that a\ < 2;=o T j < a if or 
all k, and sup^ \tj\ — C < 1, the Schwarz-Christoffel function 

is holomorphic in the upper half plane and satisfies 

| z |i-Ef=o(i-i/x>r; minlsi^ai^ + (1 - a x )p), sin(or 2 7r + (1 - a 2 )/3)} 

< 3H(\z\e ifi ) 

< fcf-ZUQ-j/XfrJ max{sin(ofi7r + (1 - a{)p), sin(a 2 7r + (1 - a 2 )p)]. 

for every f3 e [0, n], with constants depending only upon the assumed bounds and on M. 

Proof. It is well known that H(z) maps the upper half plane conformally onto the region 
bounded by a polygonal path with angles ht j at the vertices and extends continuously to 
the boundary, so that the negative real axis is mapped to itself and at w e (M k ,M k+l ) c R 
the direction of the boundary path is n £q t ;> which by hypothesis is in [a\n, a 2 n]. We will 
need to know an estimate for the direction of the image of a ray from at angle /3 e (0, n). 
At a point w on this ray, the angles fij from w to M' e R form a strictly increasing sequence 
beginning at /3o = —n + /3 and converging to 0. The direction tangent to the image curve at 
H(w) is then + 2/~T;)jfy. Now 

* k k 

P + = P + -Pi Yj T J + Z T ^ ~ 71 ~ M 

7=0 7=0 7=0 

and/3- n—fij < for all /*, so the second sum has the opposite sign to Yij T j > 0- It follows 
that 

k k 

P + Yj(~Tj)Pj <P + (*-P)J]Tj<p + (7r- P)a 2 = tt - (1 - a 2 )(7T - 0). 

j=o j=o 

To see a lower bound, write 

k k k 

p + Z (_t ^ = & + (7r - ® t ° + (7T - & Z r > + Z Tj(p -x- Pi)- 

7=0 7=1 7=0 

If Hi T 7 ^ 0- tnen so is 2* _T A*> so 

s 

/3 + Z(- T >)£> * ^ + (* - > + - i • 

7=0 

However, at any k such that Y,\ Tj < we have Y, k j=o T j(P ~ n ~ Pi> > 0, so that the same 
lower bound 

k k 

p + Z (_r ^ - ^ + (?r ~ # Z t j- 13 + (n - ® ai = n - (1 ~ ai)(n - ft 

7=0 7=0 

holds. Thus the direction of the image of the ray from at angle /3 is between n— (1 — a\){n— 
p) wdn-{\—a 2 ){n—p. In particular we may relate the integral along this ray to the integral 
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of the magnitude of the integrand along the ray. Writing h(w) — w T ° YIJ=\(J ~ W M J ) Tj 
and z = \z\e'P we have 

Hzl, 



min{sin(7r-(l-ff 1 )(7r- y 6)),sin(7r-(l-a 2 )(7r -/?))} J |/i(fe^)|<fr 
r \z\ 

< 3 h(te if3 )d(te iis ) = 3#(z) 
Jo 

(8.2) < max{sin(7r - (1 - a x )(n - (3)), sin(/r - (1 - a 2 )(n -j8))} J" |/i(fe'^)| dt. 

To proceed we need an estimate on the magnitude of the integrand h{w). For fixed z let 
£ be such that M k ~ l/1 < \z\ < M k+1/Z and write 



\ m \ = w -oh _ iL|- T * ni^p fth - n h _ jlp 

I v 7 I 1 1 I M*l A llAf/l 1 II w I 111 M'l 



= lwr To i - 



7=1 



7=1 



k-l _ t k-l 

ni^/i Tj «p(z-r,-iog|i 

7=1 7=1 



Now it is easy to see that 



k-l 



7=1 



7=*+l 



T/log 


1 - 






w 1 


T/log 


1 - 


W 1 

Mil 



7=1 



C 



J] -T.IOJ 

M+l 



w I 



j=k+i 



so that 



k-l 



k-\ 



u i — - iwr T ° — ^ i — r \M (j - k)T 

v y| I M k \ 1 HmjI I M* I 11 



MJ\ 



M k I 



with constants independent of A: (here we used boundedness of £q Tj). We may use this to 
estimate the integral in (8.2) from t = A k ~ 1/2 to t - A k+1/2 . The fact that \r k \ < 1 ensures 
integrability of the first term, and that the integration introduces a bounded multiple of M 
into the product. Letting K be the integer part of log \z\/ log M we find that 

M K k-l K-l 

\KteH dt « Y F[ M ( -'~ k)T ' > M K r~f M { i~ K)Tj . 
However the ratio of the (A + l)-th term to the £-th term is 



\ ( k-i 



M k+l Y] M (j - k - 1)T > M~ k Y] M (k - j>T ' 

7=0 J I 7=0 

by hypothesis. Thus we have the upper bound 



K k-l 



( K-l 



^ K 



K-l 



J]M k Y] M { '- k)T ' < M K \\ M (j ~ K>Tj Yj M- k<l - a ^ <M K Y] M (j - K)T i 



k=0 >=0 



V 7=0 



) k=0 



7=0 



and finally may conclude (again using boundedness of £q Tj) 

Hz|, 



[ Z \h{teH 
Jm° 1 1 



dt a M i-LUa-i/^i. 
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from which the result follows by (8.2). □ 

The construction in the preceding lemma is applicable under the hypotheses of the the- 
orem, as shown in the following. 

Lemma 8.5. Suppose f satisfies the hypotheses of Theorem 8.3. For any < e < 
| min{(/3i), (1 - f} 2 ), (J3 2 — P\)} there is a constant C3 = c^(e) and a function F(z), holo- 
morphic on the upper half plane, such that for all f3 € [0, n] 

c 3 /(k|) sin(08i + e)(;r -/?)) < %F(\z\e ifi ) < f(\z\) sm((/3 2 - e)(n-0)). 

Proof. By hypothesis we have 

ci log/(M')-log/(M'-') c 2 
Pi ~~ 1 77 - 1 77 - "2 + 



log M log M log M 

Choose M sufficiently large that both and are smaller than e. Then 

log f(M') - log f{M'~ l ) 

(8.3) 0<Bi-e<— ^- — - ^ -<B 2 + e<l. 

logM 

By making M larger if necessary we also may require f}\ - e < ^Ipjr < fii + £, because 
log(/(l))/logM -> as M -> 00. 
Now let 

, _ log/(M> +1 )_ . >Q 



and 



; O'+DlogM' 

l-l j0 if J = 

2/o-2/i ifj=l 

tj = j(2h , - lj - l h2 ) + l h2 - I j if j > 2 



We claim that 



k 

(8.4) 1-^(1- j/k)T J = l k - l fork > 1. 

7=0 

This is easily verified when k — 1,2. For k > 2 it follows by induction, because if it is true 
up to k then the the following are true, and we must prove X = (k + 1)1^. 



k - ^(k - f)Tj = kl k -\ 
j=o 

k-\ 

k-\-YjJc-l-])Tj = {k-l)l k -2 
j=0 

k + \ 

k+l- J](k + 1 - j)Tj = X. 



If we subtract the latter two equations from twice the first we obtain 

T k = 2^4-i - (k - l)l k - 2 - X 
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from which we have the desired equality: 

X = 2kl k - 1 -(k-l)l k - 2 -T k 

= 2kl k -\ -(k- 1)4_ 2 - k{2l k -\ -l k - l k - 2 ) - h-2 + h 
= (k+l)l k . 

Our goal is to apply Lemma 8.4, for which we need estimates on £ Tj. Define a\ = 
1 -/3 2 — e and «2 = 1 —p + it is easy to check that < a\ < a2 < 1 . For ; = 0, to = 1 - Iq 
is in [ori, by the estimate following (8.3). For ; > 1 we may use (8.4) to obtain 

i i+ 1 i 

Z Tj = Z (i + 1 ~ j)r j - Z (i ~ j)T j 

j=0 j=o o 
= 1 - (i + 1)/; + ?'/;_] 

_ 1 log/(M f+1 )-log/(M') 
log M 

€ [ai,a 2 ] 

directly from (8.3). We also need that sup^ \tj\ < 1. As already noted, to 6 [aijOiL For 
j — 1 we write 

, n , M , log/(M 2 )-log/(M) 

ti = 'o - U«i - <o) = h ; tz 

logM 

which is a difference of two values from - e,p + e], so has magnitude bounded by 
2e + p — Pi ■ Similarly, for j > 2 rewrite tj as 

nj \ /, . nJ , \ \ogf{Mi)-\ogf{Mi- x ) log/(M^')-log/(M^) 

which is a difference of the same type. Since we chose that 3e < p- P\ we conclude that 
all |Tj-|, j > 1 are bounded by p ~ Pi < F 

We have verified that Lemma 8.4 may be applied to the sequence t, to produce a 
function H(z). Since a-2 > u\ = 1 - /?2 - e > we may find a constant c' such that 
sin(j6 + ff2(7r - Pi) < c' sin(jS + a\(n — pi) for all /3 e [0, 7r]; inserting this into the estimate 
of 3H(z) in Lemma 8.4 yields 

| z |i-Z*o(W/K)r/ sin ^ + a^^p)) < SHiW?) < \ z \ l -ZU l -j' K > T i s in(yS + ai(n~P) 
with the integer part of log \z\/ log M. However we determined in (8.4) that 

i Z |i-zjLo(w/«T./ _ | z |fe-i _ | z | tziZsM - (f(M K )y 



Moreover (Ar+ 'ffi~g M £ [1-^T, 1] and f{M K ) lKK+l) is bounded, as is the ratio of f(M K )/f(\z\), 
so we have in fact 

f(\z\) sin(fi + a 2 {n - P) < 3//(|z|^) < f(\z\) sin(/? + a x {n - p)) 

which may be rewritten as 

f{\z\)sm{{\-a 2 ){n-P) < 3HQz\e^) <; /(|z|) sin((l - ai )(n -ft)). 

Hence there is some constant c such that F(z) = -iH(z)/c satisfies the conclusion of the 
lemma. □ 



RESOLVENT ESTIMATES AND FRACTAL BLOWUPS 



35 



9. Estimates away from the negative real axis. 

In previous sections we have estimated piecewise eigenfunctions and the resolvent for 
real positive values of A. In this section we combine them with a weak decay estimate 
for the resolvent and use the Phragmen-Lindelof theorems to obtain decay estimates in a 
sector in C. This gives bounds on the resolvent everywhere away from the negative real 
axis. 

Our weak decay estimate for the resolvent relies on knowing L°° bounds for projection 
kernels onto eigenspaces. To obtain them we first prove a modified version of Theorem 
4.5.4 of [9], which relates L°° and L 2 norms for eigenfunctions. 

Theorem 9.1. If u is in the span of the eigenfunctions with eigenvalues not exceeding A 
then \\u\U < A S/2(S+1) |M| 2 . 

Proof. Let u = £ a j<Pj( x )> where the <f>j are an orthonormal set of eigenfunctions with 
eigenvalues -Aj, and < Aj < A for all j. The main point of our argument is that u does 
not vary very much on any cell F W (X) for which r w fi w A < i. Specifically, if F W (X) is a 
cell and h w is the harmonic function on X with h w {p) — u(F w (p)) for all p e Vq, then the 
difference between u and h H . on F W {X) can be obtained by integrating against the Green 
kernel G (0, (jc,y). Since we know ||G®(jc,3')||co < 1 we have 

\\u o F w - K ° Fw|L = |JV°W)(-A(« o FJfy)) dpty)^ 

< \\a(u o fjHj 

= r w /z w ||(Aw) ° ^w||i 
= r 4 Au L ( FAX)) 

(9.1) ^/1HI M ))" 
If is a partition of X into cells with r w fi w A < \ then 

(9.2) maxr w jiJ|Aw|| = maxr w p w (V a 2 /l 2 ) 1 < maxr w p w A||M|| 2 < h\u\\ 2 . 

Now if we make the above decomposition on each cell of we obtain a piecewise harmonic 
function h by setting h o F w = h H . o F w and by (9. 1) and (9.2) 

(9.3) \\u - Alio, < maxr w /iy 2 ||AM|| L2(F ^ (x)) < maxr WJ ui, /2 ||AM|| 2 < ^ max//,;, 1/2 ||M||2 
However (9.1) and (9.2) also yield 

II" - ,l Wl = " m ~ /2 Hi 2 (f„(X)) = 2^"" ° Fw ~~ hw ° Fw ^l 2 

W W 

< 2^/x w \\u o F w - h w o F w 

W 



I.Y] 

2 



|2 



(9.4) <max^.||A«|f < j|M| 2 . 
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Moreover the L°° and L 2 norms of h w o F w are comparable because the space of harmonic 
functions is finite dimensional. As in Lemma 4.5.5 of [9] we conclude 

HAIL < Z H /! »' ° F "ll~ 5 Z m » ° Fwlb = Za 1/2||/iJl2 ^» = ™x/^ 1/2 Plb 

www 
3—1/2 

and by (9.4), p|L < | max, v // w . ||m||2- Combining this with (9.3) gives 
IML < IWL + imax^ 1/2 ||«|| 2 < 2max^, 1/2 || M || 2 . 

I w 

We required only r w /x w A < i, so may choose it such that /4f +1 ^ s = r w /i w > A -1 , from 
which ^, 1/2 <A S/2(S+1) . □ 

Theorem 9.2. If{cf>j} is a finite orthonormal set of eigenfunctions with eigenvalues not 
exceeding A and {bj} is a set of complex numbers with all \bj\ < B then 



IE^HLw BAS/(S+1) 



(XXX) 

J 

Proof. Fix y e X and take the L°° norm with respect to x. We see that <pj(x)(pj(y) is in 
the span of the eigenfunctions with eigenvalues at most A, so by Theorem 9. 1 it suffices to 
compute the Lr norm with respect to x. Since the (pj(x) are orthonormal and appear with 
coefficients <f>j(y), we have 

a~ s/(S+ 1I2>«HL w *2>^ 2 * IE w 2 IL * lE^HUx*) 

i j j 

Taking the supremum over y gives 

IE^HLw aS/(S+1) 

because the functions are continuous and X is compact. However the the fact that the <pj 
are real-valued then ensures 

a/2 /v -, , -1/2 



|2 * B Z i^-w^wi * B (Z M*) 2 ) (Z w 2 ) ~ baS/(SH 



f+i) 

; j ; j 



□ 



Remark 3. The estimates in Theorems 9.1 and 9.2 are sharp in the sense that there are 
fractals on which they are achieved. Indeed, on certain fractals with sufficient symmetry 
it is also known that there are individual eigenfunctions of eigenvalue A with support in 
cells of size A~ s t (s+l) when A is large. For these to have L 2 norm 1 they must then have L°° 
size A S/2(S+1 \ and the product (f> j(x)(f> j(y) then has size A S/(S+1) (see Theorem 4.5.4 of [9]). 
We do not know whether they are sharp on all pcfss sets, though it is not difficult to get a 
lower bound from the Weyl law in Proposition 3.1. According to that result, the number of 
eigenfunctions with eigenvalue at most A is bounded below by a multiple of A s ^ s+1) if A is 
sufficiently large. Calling them <pj we see <pj(x)4>j(y) has L 2 (XxX) norm bounded below 
by a multiple of A S/2(S+1) , and hence U°{X x X) norm that is at least this great. Hence as 
a converse to Theorem 9.2 we have the existence of a set of eigenfunctions satisfying the 
assumptions and for which there is x at which £ ; (f>j(x)<pj{x) > \ s / 2 ( s+l \ The consequence 
for Theorem 9. 1 is that setting aj = (f)j(x) we have the L°° norm of u = Yij a j<Pj( x ) is equal 
2jfl 2 <; A s ^ s+l \ while its L 2 norm is (Zy'tfp ' 2 > so u satisfies the hypotheses and has 

IML > a s ^ s+1) ||m|| 2 . 
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Lemma 9.3. For z € A a , the piecewise eigenfunction rfp satisfies 



K :, W|| L „ (Z) <(l+tan(f)). 



2' 

Proof. Recall ( p is the harmonic function on X with value 1 at p e Vo and at the other 
points of Vo. Let {(f) j] be an orthonormal eigenfunction basis of L? with corresponding 
eigenvalues {—Aj} repeated according to multiplicity. It is easy to verify (see Lemma 3.3 
of [6]) that if the L? expansion of £ p in this basis is £ p — Y,j bj4>j then we have g p - rfp = 
Z 2,i "^T^h w i trl ^ 2 convergence of the function and its Laplacian. We also have the L 2 
expansion rff = j^'f'j- Now we compute at z e A a 

(9.5) r,f( x) - rfpM - (|z| - z) £ + ^ + A 

j 

We break up the sum into pieces of the form Aj e [2 k , 2 k+l ) and estimate each piece using 
Theorem 9. 1 to obtain 

i^(x) - ,r «iu * n - 4 z w fe+ Aj£ + ^ L 

<li z i_ z |y 2 (™ 2 ' s+i )|| y — x -h — 



|£ 2 (* + 1)S/2(S + 1)/ ^ I^A'I 2 \' /2 

I^AI 2 



yUs^ max Ezsy y W \ 



Now for k > Icq — log 2 \z\ + 1 we have Aj > 2\z\ so \z + /l/| > 4- > 2* _1 and the maximum in 
the above is at most 8|z|2~*. For Htowe instead bound the maximum using the fact that 
z £ A a . Clearly we have |z + A\ > \z\ sin a, but this is a poor estimate when a < |, where 
\z + A\ > \z\. Instead observe that sin a > cos | for a e [f ,tt] and so |z + A/I > |z| cos | for 
all e [0,7r]. Also ||z| — z| < 2|z| sin |, so the maximum is bounded by 2tan |. Inserting 
these bounds for the specified k and using the Cauchy-Schwarz inequality gives 

\\nf{x) - r]f\x)\U 

<2tan(^)JV S/2(S+I) ( J] i^ A Jl) 1/2 + 8|z|^2^ +1 )2W £ 

2 i<io 2*</L-<2* +1 |kl + &]\ k>ko 2 i <,l / <2 1 '+ 1 |kl + Aj| 

<2tan(f)(^2^-") I/2 ( £ ^ 2 f + 8|z|g 2^2^ f\ £ 



< (2tan(^)2 w/2(s+1) + 8|z|2™ +1 >2-*")(£ ^ ) 
<(4tan(^)|z| s/2(s+1) + 16|z| s/2(s+1) )||^ 



(1:1) 1 1 



where in the final step we used that fco = log 2 |z| + 1 and the known L? expansion of 
T]p Z<l \ The latter can be estimated using (5.1) from Theorem 5.1. This shows rf^ has 
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exponential decay with scale k(\z\), so that its L? norm is dominated by integrating over the 
cell containing p and of measure |z| _ ' s '(' s+1 \ Since the function is bounded by 1 we find 
H'^IL ~ |7p s / 2 (' s+1) . Again using (5.1), this time as an L°° bound on 77^", completes the 
proof. □ 

Lemma 9.4. For z € A a , the Green kernel for a pcfss fractal with regular harmonic struc- 
ture satisfies 

ll^y)|U x ,^(l + kir 1/(5+1) (l + tan(f)). 

Proof. Expanding the resolvent kernel at both zeA„ and at \z\ e R with respect to the 
basis \<pj) of eigenfunctions from the previous theorem we find 

g^W) - G^ y) - £(-J_ - ^M)Hy) - Z (z + ^ +A]) H^- 

As before we break up the sum, taking one piece where Aj < 4\z\, and the rest to be 
of the form Aj e [2 k , 2 k+l ), k > ko — log 2 |z| + 1. For each of the the latter we have 
|z| + Aj > \z + Aj\ > if > 2* _1 so the coefficients multiplying (f>j(x)cpj(y) are bounded by 
8|z|2" 2 \ Applying Theorem 9.2 we find 

<m> II S (J^/ ^HL^ 2 ^''- 1 '"" 1 ' 

2*|z|</l^2* +I |z| v JJ ' 1 

For the piece where Aj < 4\z\ we use \z + Aj\ > \z\ cos ^ and ||z| - z\ < 2\z\ sin f for all 
a 6 [0,7r], because z e A a . The coefficients are then bounded by 2|z|~' tan |, and applying 
Theorem 9.2 gives 

II Z } TT^TT^ Mtoi ^2| Z r 1 ta n (?)kl S/(S+1) = 2kr 1/(S+1) tan(^). 
1 (Z + ^/Xlzl + Aj) Hoc 2 2 

Combining this with (9.6) we have 

< 2|zr 1 ^ +1 » tan(-) + 8|z| Z 2" 2i 2* +1 ^^ + " 

< |zr 1/(s+1) tan(|) + 8 \ z \2- 2k °2^ +l)s/( < s+1) 



">(l+tan(|)) 



-l/(S + l)f 

2' 

However the bound can be improved for small z, because if 2\z\ < Aq, the smallest Dirichlet 
eigenvalue, then the first sum is empty and the second begins at log Aq rather than log \z\, 
giving abound independent of z- We may therefore replace |z| _1 ^ s+1) with (1 + |z|) _1 ^ s+1) in 
the estimate. The result now follows by using (6. 1) from Theorem 6. 1 to bound G ( '*(x,y). 

□ 

Remark 4. A small modification of the method in the previous proof provides bounds off 
suitably small neighborhoods of the eigenvalues. If Dj is the disc radius 8Aj ^ centered 

at -Aj then outside UjDj we have \z + Aj\ > S^A 1 / 1 * and using this for those Aj < 4\z\ 
rather than the estimate in the above proof we see that 

\\G®(x,y)\\ sa-V-WDoff UjDj 
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Moreover the Weyl estimate in Proposition 3.1 guarantees that the number of eigenvalues 
of size at most A is bounded by a constant multiple of A ,s/(,s+1) and therefore a suitably 
small 6 depending only on the fractal and harmonic structure ensures UjDj contains at 
most half of any interval of the form [— 2' +1 , — 2*]. It is easy to check that this estimate and 
that in the lemma are comparable at points of the form —Aj + iSA\ . 

We complete our estimates in A a with some for the normal derivatives. 

Lemma 9.5. Under the assumptions of Lemmas 9.3 and 9.4, 

(9.7) \d n rfp\q)\ S d + kl)' /(S+1) (l + tan(|)) 

(9.8) \\d' n G^(p, X )\\ L ^ X) <(l+t^)) 

(9.9) |«G<*W| < (1 + kl) 1/(S+1) (l + tan(|)). 

Proof. Applying the Laplacian to 

G (z \x,y)-G (lzl \x,y) = Y — 0,(x)0,-(V). 

gives an L 2 convergent series. We may therefore expand in the Gauss-Green formula 

d'^G (z) (x,p)-d' i ;G m (x,p) 

= d',;(G {z \x, q) - G m {x, q)%{q) ~ (G (z) (x, q) - G (l «(x, q))d n £ p (q) 

qeV 

= j A(G®(x,y)-GM\x,yM P (y)dfi(y) 

:\z\-z)Aj 



as we saw in (9.5). The justification for exchanging the integral and sum is that \\<pj\\oa < 
A s/2(s+i) by 9 { and ^ 

is in L 2 , so the series form of the integrand is L l convergent 

uniformly in x. The same is true of the series expansion of rfp(x) - r]p(x) by Lemma 9.3, 
and combining the estimate of that lemma with (6.2) from Theorem 6. 1 gives (9.8). 
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Similarly, the fact that applying the Laplacian to rfp(x) - rj p (x) gives an L convergent 
series allows us to compute 



qeVo 

= j h{rjf{x)-r 1 f\x))i p {x)dp.{x) 



(z + Aj)(\z\ + Ay) 



<f>j{x)( p {x)d[i{x) 



<(lz|-z)£ 



( Z + ^)(|z| + /lj) 



-^(kl + ^) 2 
<2|z|||<'f £|z| 1/(S+1) 

where we used the decay estimate (5.1) as we did in Lemma 9.3. This and the normal 
derivative estimates from Theorem 5. 1 give (9.7) and together with (6.3) from Theorem 6. 1 
establishes (9.9). □ 

With the preceding weak decay estimates in hand we may prove the main result of this 
section, which gives bounds on both the functions rfp and the resolvent kernel G^ z \x,y), 
as well as their normal derivatives, at points z e C that are not on the negative real axis. 

Theorem 9.6. Suppose X is a post-critically finite fractal with regular harmonic structure. 
Then there is a constant k$ > depending only on the fractal and harmonic structure such 
that for z — \z\e'@, \j3\ < n we have 

\r]f(x)\iZ®(J3,\z\,P,x) 

\d n nf(q)\^k + i\ ms+1) ^,\z\,p,q) 

\G (z \x,y)\ < \z+ ir 1/(5+1) <B08, \z\,x,y) 
\d'„G <z \p,y)\<<S>(fi, \z\, P ,y) 
\d:;d'„G iz) ( P ,q)\ < \z+l\ ms+1) ^,\z\,P,q) 

where 

a,08Jz|,x > y) = (l + tan(^))exp(-^ l )fey)sin(^M)) 

Proof. We give only the proof for G^(x,y) because the others are entirely analogous. Fix 
x and y in X. Recall from (6. 1) that the resolvent satisfies 

(9.10) G w (x,y) < (l+Ar 1/(s+1) exp(-K 2 d kw (x,y)). 

for A e (0, oo). 

Let f(A) = K2dk(x«i')(x,y). Recalling Definition 5. 12 we may substitute into Lemma 4.2 
with k' = - log M and k = jkf to find that there are constants c\ and c 2 , depending only on 
the harmonic structure, for which 

_a_ ^( jW 0+D"/")(^,3 ; ) , » 

Cl M-' s +» < — j - < c 2 M^ 

dk(Mi-i'){x, y) 
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provided M is large enough. (Note that the latter restriction is so k{M ( j +l) "l n ) is given by 
the formula in Definition 5.12 and is not zero.) Thus (8.1) holds for f(A) with < /?i = 
jrfjjjTj < Pi — ^1 < 1 ■ Consider the product 

( Z+ 1)1/(S + D 

s(z) = — G {z) (x v) 

gK> C(l+tanf) V,y) 

on the sector A+. Lemma 9.4 implies \g\ is bounded by 1 on AJ if C is large enough, 
while (9.10) implies \g(z)\ < exp(-/(|z|^ Q; )) on the positive real axis. Applying Theo- 
rem 8.3 and multiplying out to retrieve G (c) we have for /3 e [0, a] that 

|G flde ^(x,y)| < \z + ir 1/(s+1) (l +tan(^))exp(-cK 2 d k(U} (x,y)smc(l - ^ 

In order to obtain an estimate on a general ray z = \z\e'P, ft + n it then suffices to take 
a = j(n + P), so 

\G^\x,y)\ < \z+ + tan(^))e X p(- C ^ 0z|) (x,y)sin(^^)) 

and by symmetry corresponding estimates are valid on the sector A~ obtained by reflection 
in the real axis, simply by replacing ft with \/3\. □ 

Note that in the case of an affine nested fractal we could have proved this result using 
the classical Phragmen-Lindelof theorem, because d^ix^) — (1 + X) y ^ s+ ^R{x,y) y by 
Proposition 4.3 and Definition 5.12. In this case one also obtains better constants. 

Theorem 9.7. The sequence G^ n {x,y) defined in Theorem 7.3 for a blow-up Q. ofX con- 
verges uniformly on compact subsets of Q X £1 X (C \ (—00, 0]) to the Laplacian resolvent 

®(J3,\z\,x,y) = (l +tan(^))exp(-^(|z| ) (x,y)sin(^ 7 ^)). 

Proof. Using the definition of G^l(x,y) and the bound from Theorem 9.6 we have 
\G^(x,y)\ 

< ^Jr^^z + ir 1/(S+1) OC8, r^ lM w„\z\, ^„ ° ■ ■ ■ o F ai x, F a „ o - • o F^y) 

= |z + VUVlJ 1/<S + 1)( I ) 03, rfjj^^Jzl, /v„ o ■ ■ ■ o F^X, F^ o ■ ■ ■ o F^y) 

but we recall from (7.6) that 

d k { r -^-^A){F<T„ ■ ' ' F^x, F lTn o ■ ■ ■ o F ai y) * d m (x,y) 
from which there is K(, > such that 

|G^(*,y)| < [z + rtai^Eoi. |" 1/CS+1) *0S, |zl. x. y) 

independent of n. Thus the sequence is uniformly bounded on any compact set of the type 
in the hypotheses, and since it is analytic in z it is a normal family by Montel's theorem. 
The function G ( ^\x,y) is the unique limit point because the sequence converges on the 
positive real axis by Theorem 7.3. □ 



G„ \x,y), and there is > such that 
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10. Estimates for other kernels 

One of the main purposes for proving estimates on the resolvent kernel is that we can 
obtain estimates of the kernels of other operators from it using functional calculus. Specif- 
ically, if r is a contour in C that surrounds the spectrum {-Aj} of A and h{z) is analytic in 
a neighborhood of F and its interior then we define 



( 10. 1) H(x, y) = ^-f T G iz \x, y)h(z) dz 

provided h is such that the integral converges. Under reasonable assumptions this should 
be the kernel of h(A). At times it is useful to think of the kernel of h(A) as H{x,y) - 
h(-Aj)(f>j(x)<pj(y), so we begin by describing sufficient conditions for these to coincide. 
It will be convenient for us to work with contours that lie in a sector A a as in Section 8. 

Lemma 10.1. Suppose T is a contour in A a that surrounds the spectrum of A, that h(z) is 
analytic in a neighborhood ofT and its interior, and that there is some a > 2 (s+i) suc ^ mat 

(10.2) I \h(z)\\zr l \dz\ < oo. 



I' 



If in addition J^j \h(— Af)\ 2 < oo then both H(x,y) and H(x,y) converge in L 2 (X X X) to a 
common limit. If the stronger condition 

(10.3) y 2 *s/(s+i) sup \h{-Aj)\ < oo 

£ 2 k <Xj<2 k * 1 

holds then both H(x,y) and H(x,y) converge in L°°(X XX) to a common limit. 
Proof. The Cauchy formula guarantees that 

H(x,y) = ^ h(-Aj)<Pj(x)<Pj(y) = ^ ^ f dz<f>j(x)<Pj(y) 
j j r ~ j 

so we must prove that we can exchange sum and integral to obtain 

f Y^-<t>j(x)<Pj(y)dz = H(x,y) 
2m J T t-f z + Aj 

j 

where the convergence is in L?(X x X) for p — 2 or p = oo. So it suffices to prove 
convergence of 

(10.4) [J] h(-Aj)<f>j(x)<f>j(y)\\ Di{xxX) and, 

j 

(10 - 5) llXS^^wHLcM,- 

For (10.4) we have 

Ehi-Aj^jix^jfyi = (Y)K-Aj)\ 2 ) m 



MZ_J IIL (XXX) ^ 2 i <,l J <2'+ 1 

where the latter is from Theorem 9.2. Summability of these are the conditions given in the 
hypotheses for p — 2 and p — oo respectively. 
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For (10.5) we use the integral Minkowski inequality and Theorem 9.2 to obtain 

llf yJ^L Hx)4>j{y)dz \ 



< 



f \Kz)\Y 2kSKS+l) su p lz + ^l -1 kfel 

Jr t 2 i <i i <2 1+1 



1/2 



P = 



But now |z + A/| 1 < tan(|) min{|z| .A- 1 }. For a > and using the Weyl law of 

Proposition 3.1 

G \Z + ^r 2 )' /2 £ feT 1 tan(f )E 2^/^'>2^) 1/2 < c^f" 1 tan(f ) 

and also 

£ 2^< s+1) sup k + Xf l < \zr l tan(f )(2 2^ <s+1 >2- 2fli ') 1/2 < cjzl " 1 tan(|) 
so that we have for both /? = 2 and p = oo and any a > there is c fl so that 

II f X-^tW^)^! ~ c « f woikr 1 1*1 

IIJ r Z T J z + / l. Nl"(ZxX) J r 

and finiteness of this integral is enough to ensure (10.5). □ 

We can then obtain L°° bounds for the kernel H(x,y) by inserting our estimates on 
G (z) (x, y) into (10.1) and integrating over F. How best to do this depends very much on the 
function h(z), so rather than attempt to formulate a general theorem we simply give some 
examples of how this can be used. There are many other examples that could be treated in 
a similar manner to those given below, such as exponentials of complex powers of A, some 
of which we may return to in future work. 

Example 1 (Heat Kernel). The heat kernel is defined for t > by 

Pr (x,y) = Y j e- A '<4,j(x)cf>j(y) 
j 

so in the above notation it is H(x,y) for h(z) = e' z . To estimate it we introduce a suitable 
family of contours. 

Fix t > and a e (0, n), and let Y a t be as in Figure 1. It consists of the arc of the circle 
t\z\ — 1 in the sector A a , as well as the rays z = k|e ±!tv , |z| 6 [r 1 , oo), traversed so as to wind 
once around any point on the negative real axis. 

Observe that on the circular arc of T^j^t we have \e ,z \ < e, while on the rays \e' z \ < e~'l ; l/ 2 
where \z\ > f _1 . It is then easy to verify (10.2), for example with a = 1. We also note that 
h(-Aj) = e~ a ', so that (10.3) is just finiteness of £ t 2 ks ^ s+ ^e~ ak . We conclude that 

Pt(x,y)=^—, \ G (z \x,y)e zt dz 
and can apply our estimates of G (z> to estimate p t . 
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Figure 1 . The contour T aJ 



Theorem 10.2. On a pcfss set with regular harmonic structure there is k~i > such that 

\p,(x,y)\ < r s/(S+1) exp(-K 7 4(ri)(*,:y))- 
Proof. We make the trivial computation 

f G (z) (x,y)e zt dz < sup \G (z) (x,y)\ \ \e z '\ \dz\ < - sup \G (z \x,y)\ 
Jr„., zer aJ Jr,,,, f zer„_, 

and use the bound \G (z) (x,y)\ < \z + |z|, JC.y) in the sector A ff/4 from Theo- 

rem 9.6, and the fact that \z\ > r 1 on the contour. □ 

This is the same upper bound as in Theorem 1.1 of [4], though their result is stronger 
because it also includes a lower estimate from which it follows that this upper bound is 
sharp. 

Note that in the case of an affine nested fractal, where by Proposition 4.3 and Defini- 
tion 5. 12 we have d m (1 + Ay KS+1) R(x,y) y this becomes 

\p t (x,y)\ < t- SKS+1) exp(-K 7 R(x, y y(l + rV /(S+1) ) 

whichshould be compared to the upper bound in Theorem 1.1 of [3] orTheorem6.1 of [4]. 
Our y corresponds to if^y in the latter reference. 

It is also worth mentioning that the above method gives estimates on the derivatives of 
the heat kernel, simply because differentiation in t introduces a polynomial power of z into 
the integrand. 

Example 2 (e wA ,w € C). Consider h w (A) for h w (z) = e wz and w e C. In order that 
that (10.3) holds it is necessary and sufficient that w = \w\e i/3 with \{3\ < f . For (10.2) 
we must be careful to select a suitable contour. If we use T a j then on the radial parts of 
the contour we have \h(z)\ = e^ e ' with the ± sign determining whether we are on the 
upper or lower ray. Then (10.2) holds for some (indeed any) a > ^ iff ft + a e (|,tt), 
meaning that rotating the boundary rays of the sector by /3 still leaves them in the left half 
plane. This is true iff l^-ffl < /3. With these constraints, following the proof of the previous 
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theorem with t — \w\ we have a bound of the form 

\H w (x,y)\ < M- s/(s+1) (cot |0) exp(-K 5 4(M-i)(*,)0 sin(e|0)). 

Example 3 (Complex powers). If we take the usual logarithm with branch cut on the neg- 
ative real axis then h(z) = exp(wlog(-z)) = {—z) w is analytic in a neighborhood of the 
spectrum of A. We cannot use one of the T a j contours, but can take the rays at angle ±a 
from radial distance ol A ° , , where -/In is the smallest Dirichlet eigenvalue, and connect 

2| COSQ'I u ° ' 

these with a vertical line through -| on the real axis. This change does not affect the 
convergence of (10.2), which depends only on the behavior of h(z) on the rays. On these 
\h(z)\ is dominated by e li™(«OI(^-l«l)| z |Re(..0 ; so (10 2 ) i s satisfied if a + Re(w) < for some 
a > 2 (s+i) ' nence we need only Re(vv) < - 2( ^ +1) . Applying the same estimate to h{—Aj) 
we see from the Weyl estimate that 

^ \h(-Aj)\ 2 < C J] ^; Re(H,) < C ^ 2^ (s+1 >2 2Re <^ +1 > < oo 

j j k 

when Re(w) < - ^J +Y . ■ This is sufficient to conclude that H(x,y) and H(x,y) converge to 
the same limit in L 2 (X x X). If we want convergence of H(x,y) also in L°°{X x X) then we 
instead need Re(vv) < — j^j ■ 

In estimating \H(x, y)\ along the rays using Theorem 9.6 we do not need any assumptions 
on w. We have 

f \h(z)G (z \x,y)\\dz\ 

jRay 

< ^Mj, + tan(^)) J k |^)-i/( S+ D exp(-^^ (|z|) (x, y) sin(c^)) \dz\ 

however dkq z \)(x, y) > (1 + \z\) l ^ <s+l> R{x,y) from the discussion preceding Lemma 4.2 and 
Definition 5.12, so this integral is bounded by a Gamma function. No matter how large 
a negative value of Re(vv) we have, we cannot obtain anything better than \H(x,y)\ ^ 1 
because the contour integral for H(x, y) includes integration along the vertical line on which 
the most we can say is that h(z)G (z) is bounded by a constant depending on the fractal and 
harmonic structure and the length is bounded in the same way. Abusing notation to use 
F( ) for the Gamma function as well as the contour we find that our bound is 

|BM sl+e ,«*»( lM £^ ) )( csc(e ^ ) ^) s *-'-« r(s+(s+1)Re( 

If | Im(w)| is larger than S + {S + 1) Re(vv) then taking a so | Im(w)|(7r - \a\) = S + (S + 
1) Re(vv) we find that the terms involving a cancel with the Gamma factor to leave only 
| Im(w)| (s+1)(Re(w)+1) . If Im(w) is smaller than this then we may take a = n - 1 and bound 
by a multiple of e lm(w) r(s + (S + l)Re(w)). In either case we get no worse than 

\H(x,y)\ <1+(S+(S + D^D^+lX^wHl^yji-^+iXReCwHl), 

We further remark that if the fractal is affine nested, so that dk(x, y) e ky R(x, y) y for 
some y e \ j] as in Proposition 4.3, then we can improve the bound to 

|H(jt,y)| < 1 +(S +(S + DlwD^ + iKReW+iHr-D/y^^^d-C^iXReOvHlw/y^ 

Some results related to these appear in [7], where they were obtained for affine nested 
fractals by recognizing the kernels as being of Calderon-Zygmund type. 
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